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ABSTRACT 


THE  APPLICATION  OF  A  DETERMINISTIC 
SPECTRAL  DOMAIN  METHOD  TO  THE 
ANALYSIS  OF  PLANAR  CIRCUIT 
DISCONTINUITIES  ON  OPEN 
SUBSTRATES 


by 

JAMES  STUART  MCLEAN,  B.S.,  M.S. 

SUPERVISING  PROFESSOR:  Tatsuo  Itoh 


A  deterministic  formulation  of  the  method  of  moments  carried  out  in  the 
spectral  domain  is  extended  to  include  the  effects  of  two-dimensional  two- 
component  current  flow  in  planar  transmission  line  discontinuities  on  open 
substrates.  The  method  includes  the  effects  of  space-wave  and  surface-wave 
radiation  through  the  use  of  the  exact  spectral  domain  Green’s  function.  The 
procedure  and  formulation  of  the  method  are  described  in  detail.  Also,  tech¬ 
niques  used  to  increase  the  numerical  efficiency  are  described  in  detail.  The 


in 


method  is  used  to  determine  accurate  circuit  models  of  three  types  of  planar 
circuit  discontinuities  on  open  substrates:  microstrip  open-end  discontinuities, 
slotline  short-circuit  discontinuities,  and  microstrip  gap  discontinuities.  The 
analysis  is  then  applied  to  gap-coupled  resonators.  The  coupling  between  cas¬ 
caded  gap  discontinuities  is  shown  to  be  significant  when  the  substrate  is  elec¬ 
trically  thick  since  surface  wave  excitation  is  strong.  Possibilities  for  further 
applications  of  the  method  to  more  complicated  discontinuities  are  discussed. 
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Chapter  1 


INTRODUCTION 


1.1  Motivation 

Planar  circuits  and  antennas  are  attractive  because  of  low  cost,  manu¬ 
facturability,  small  size  and  light  weight.  Recently,  integrated  planar  front-end 
systems  combining  both  antennas  and  circuits  have  been  developed.  Because 
planar  circuits  are  extremely  difficult  to  modify  and  tune,  very  accurate  mod¬ 
els  are  required  for  the  design  process  in  order  to  avoid  costly  design  itera¬ 
tions.  The  modelling  of  unshielded  open  planar  circuits  is  complicated  by  the 
possibility  of  surface  wave  and  space  wave  radiation.  In  the  case  of  open  pla¬ 
nar  transmission  line  structures,  discontinuities  can  generate  space-wave  and 
surface-wave  radiation.  This  may  be  desirable,  as  in  the  case  of  planar  anten¬ 
nas.  For  example,  in  Figure  1.1  [lj,  a  planar  antenna  array  is  shown  which 
makes  use  of  microstrip  open-end  discontinuities  as  radiating  elements.  In  the 
case  of  most  circuit  applications,  radiation  is  undesirable  because  it  can  cause 
loss  and  extraneous  coupling.  Extraneous  radiation  and  coupling  may  also  ad¬ 
versely  affect  the  performance  of  antennas.  In  the  antenna  shown  in  Figure  1.1, 
radiation  ideally  occurs  from  the  open-end  discontinuities.  However,  radiation 
from  other  parts  of  the  circuit  can  occur  as  can  coupling  between  radiating  ele¬ 
ments  through  both  space  wave  and  surface  wave  propagation.  It  is,  therefore, 
necessary  to  include  these  effects  in  the  design  of  the  antenna. 
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Figure  1.1:  5-Element  Microstrip  Combline  Travelling- Wave  Antenna 

A  grounded  dielectric  substrate  can  support  both  transverse  mag¬ 
netic  (TM)  and  transverse  electric  (TE)  surface  wave  modes  of  propagation. 
For  electrically  thin  substrates  ( d  <  .01A),  surface  waves  are  not  strongly  ex¬ 
cited  [2].  However,  in  the  case  of  electrically  thick  substrates  such  as  GaAs 
millimeter-wave  integrated  circuit  (MMIC)  substrates,  the  effects  of  surface 
waves  are  more  prominent.  The  cutoff  frequencies  of  the  even  TM  and  the  odd 
TE  surface  wave  modes  of  a  grounded  dielectric  slab  are  given  by  [3] 

r  _  nc° 

hn  ~  2dv^n:* 

where  d  is  the  thickness  of  the  substrate,  e/?  is  the  relative  dielectric  constant 
of  the  substrate,  Co  is  free  space  wave  velocity  and  n  is  the  order  of  the  surface 
wave  mode  (n  =  0, 2, 4, . . .  for  TM  and  n  =  1, 3, 5, . . .  for  TE).  In  Table  1.1, 
the  cutoff  frequencies  for  the  4  lowest  order  surface  wave  modes  are  given  for 
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mode 

cutoff  frequency 

25  mil  Alumina  (ch  =  9.9) 

25  mil  GaAs  (e/*  =  12.8) 

TM0 

DC 

DC 

TE, 

39.590  GHz 

34.380  GHz 

tm2 

79.180  GHz 

68.770  GHz 

te3 

118.77  GHz 

103.15  GHz 

Table  1.1:  Cutoff  Frequencies  for  Surface- wave  Modes  of  Conductor-backed 
Dielectric  Slab 

two  widely  used  substrates.  Because  the  TMq  surface  wave  mode  has  no  cutoff 
frequency,  there  is  always  at  least  one  surface  mode  which  can  propagate  on 
an  open  substrate.  At  millimeterwave  frequencies,  several  surface  wave  modes 
may  propagate.  For  example,  in  the  upper  Ka  band  (33  to  36  GHz),  both 
TMq  and  TE\  may  propagate  on  a  conductor  backed  25  mil  GaAs  substrate. 
Surface  wave  excitation  can  result  in  losses  and  extraneous  coupling  in  open 
planar  cicuits.  The  excess  coupling  may  or  may  not  be  desirable  depending 
on  the  application.  However,  because  these  effects  become  prominent  at  high 
frequencies,  it  is  necessary  to  develop  full-wave  models  so  that  they  may  be 
accounted  for  in  the  design  of  planar  circuits  and  antennas. 

1.2  Previous  Work 

Discontinuities  in  planar  transmission  lines  have  been  the  subject 
of  numerous  analytical  investigations.  Most  of  these  analyses  make  use  of 
quasi-static  approximations  [4]-[12].  These  approximations  limit  the  analysis 
to  frequencies  at  which  the  dimensions  of  the  circuit  are  small  compared  to  a 
wavelength,  and  the  effects  of  radiation  may  be  neglected.  This  dissertation 
is  concerned  with  the  full-wave  analysis  of  discontinuities  in  planar  circuits  on 
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open  substrates.  By  full-wave  analysis  we  mean  that  the  method  of  analysis 
makes  use  of  an  exact  Green’s  function1  for  a  current  source  above  a  grounded 
dielectric  slab.  Thus  the  analysis  includes  the  effects  of  space- wave  and  surface- 
wave  radiation.  It  does  not  include  the  effects  of  conductor  loss  or  conductor 
thickness. 

Some  work  has  already  been  done  in  the  area  of  fullwave  analysis 
of  planar  transmission  line  discontinuities.  All  of  the  work  discussed  below 
has  in  common  the  use  of  the  exact  Green’s  function  for  the  grounded  di¬ 
electric  slab.  The  work  may  be  divided  into  two  categories,  eigenvalue  and 
deterministic  approaches.  The  eigenvalue  approach  involves  determining  the 
eigenfrequencies  of  a  resonant  structure  and  then  extracting  the  parameters  of 
the  discontinuity.  For  example,  a  long  rectangular  resonator  may  be  analyzed 
in  order  to  determine  the  characteristics  of  an  open-end  discontininuity  [14]. 
First,  the  resonant  frequencies  of  the  resonator  are  determined.  Then,  from  the 
resonant  frequencies  and  the  physical  length  of  the  resonator,  the  capacitance 
and  conductance  of  the  open-circuit  discontinuity  can  be  determined.  The  de¬ 
terministic  approach  involves  making  use  of  a  source  formulation  to  excite  a 
system.  In  the  case  of  the  open-end  discontinuity,  an  incident  wave  may  be 
used  as  the  source  formulation,  and  the  complex  reflection  coefficient  of  the 
open-end  may  be  solved  for  directly.  Each  method  has  specific  advantages  and 
disadvantages  as  described  below. 


Several  different  Green’s  functions  have  been  derived.  See  [13]. 
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1.2.1  Eigenvalue  Formulation 

•  The  variational  formulation  reduces  error  introduced  by  inadequacies  in 
basis  functions. 

•  The  determinant  search  is  numerically  expensive.  In  the  case  of  a  ra¬ 
diating  structure,  the  eigenvalue  becomes  complex  and  the  search  two- 
dimensional. 

•  Network  parameters  are  not  directly  determined.  Instead,  they  must  be 
inferred  from  resonant  frequencies  and  Q  values.  This  is  inconvenient  and 
can  be  numerically  unstable. 

1.2.2  Deterministic  Formulation 

•  Network  parameters  may  be  directly  calculated.  Deterministic  formula¬ 
tions  allow  for  the  direct  calculation  of  scattering  parameters  or  short- 
circuit  parameters 

•  A  determinant  search  is  not  required.  The  deterministic  formulation 
requires  only  the  filling  and  inverting  of  a  single  matrix. 

•  The  formulation  for  the  network  parameters  is  not  variational.  Therefore, 
the  error  introduced  through  inadequate  basis  functions  is  not  reduced. 

Jackson  and  Pozar  [15]  analyzed  microstrip  open-end  and  symmetric 
gap  discontinuities  using  a  moment  method  in  the  spectral  domain.  They  ne¬ 
glected  transverse  current  flow  on  the  strip  and  mode  conversion.  They  made 
use  of  a  travelling-wave  source/load  formulation.  Katehi  and  Alexopolous  [16] 
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also  analyzed  open-end  and  symmetric  gap  discontinuities  using  a  moment 
method  in  the  space  domain.  Their  analysis  was  performed  using  a  voltage- 
gap  generator  source  formulation.  They  used  approximations  similar  to  those 
used  by  [15].  Boukamp  and  Jansen  [17]  have  analyzed  open-end  discontinu¬ 
ities  in  a  covered  but  laterally  open  environment  using  a  moment  method  in 
the  spectral  domain.  Their  analysis  included  the  effects  of  transverse  current 
flow  on  the  strip  but  not  the  effects  of  mode  conversion.  They  made  use  of  a 
travelling- wave  source  formulation  similar  to  [15].  Yang  and  Alexopolous  [18] 
analyzed  open-end  and  gap  discontinuities  in  an  open  microstrip  with  a  dielec¬ 
tric  superstrate  using  a  spectral  domain  moment  method  and  a  travelling- wave 
source  formulation.  Drissi  et  al  [19,  20]  analyzed  symmetric  gap  discontinuities 
in  open  microstrip  using  a  space  domain  moment  method  applied  to  the  mixed 
potential  integral  equation.  Although  it  is  not  made  clear  in  their  paper,  they 
seem  to  have  included  the  effects  of  two-component  current  flow  and  mode 
conversion.  They  used  a  voltage-gap  source  formulation  similar  to  [16]. 

Jackson  [21]  has  analyzed  open  circuit  discontinuities  in  coplanar 
waveguide.  Yang  [22]  has  analyzed  short-circuit  slotline  discontinuities  using 
a  moment  method  in  the  spectral  domain.  He  assumed  a  purely  transverse 
electric  field  in  the  slot. 

1.3  Organization  of  This  Report 

The  remainder  of  this  report  is  organized  as  follows.  In  Chapter  2, 
the  theory  underlying  the  method  is  presented.  The  details  of  the  numerical 
computation,  including  the  selection  of  basis  functions  and  techniques  to  speed 
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up  the  numerical  integration,  are  given  in  Chapter  3.  In  Chapter  4,  the  method 
is  applied  to  the  microstrip  open-end  discontinuity.  In  Chapter  5,  the  method 
is  applied  to  the  short  circuit  slotline  discontinuity.  In  Chapter  6,  the  method 
is  applied  to  the  microstrip  gap  discontinuity  and  the  gap-coupled  resonator. 
Finally,  in  Chapter  7,  some  conclusions  are  drawn  and  suggestions  for  improve¬ 
ments  to  the  method  are  made.  Also,  further  applications  of  the  method  are 
proposed.  The  Green’s  functions  used  in  the  analysis  are  derived  using  the 
immitance  approach  [13]  in  Appendices  1  and  2. 


Chapter  2 


THEORY 


In  this  chapter  we  present  the  theory  underlying  the  deterministic 
spectral  domain  analysis.  Only  the  theory  common  to  all  of  the  analysis  will  be 
presented  here;  the  details  of  the  analysis  as  it  applies  to  specific  discontinuities 
will  be  presented  in  subsequent  chapters.  We  begin  by  developing  the  integral 
equation  which  will  form  the  basis  for  all  of  the  analysis  in  this  study. 

2.1  Integral  Equation  Formulation 

In  this  section,  we  develop  the  integral  equation  appropriate  for  the 
analysis  of  microstrip  line  and  microstrip  discontinuities.  Referring  to  the  coor¬ 
dinate  system  in  Figure  2.1,  the  electric  field  in  the  plane  y  =  d  due  to  currents 
localized  in  the  plane  y  =  d  can  be  expressed  as  the  convolution  of  a  source 
current  and  a  Green’s  function.  The  components  of  the  Green’s  function  rep¬ 
resent  the  x  and  z  directed  electric  field  due  to  a  z-directed  and  an  x-directed 
infinitesimal  horizontal  current  element  (Hertzian  dipole)  over  a  grounded  di¬ 
electric  slab  of  infinite  extent.  Since  the  Green’s  function  is  actually  a  dyadic, 
we  write  two  equations,  one  for  the  z-directed  electric  field  and  one  for  the 
x-directed  field. 
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9 


y 


Figure  2.1:  Coordinate  System  for  Dielectric  Slab 


Ez(x,d,z)  =  J  J  Zzz(x  ~  x0,  z  -  z0)J2(x0,  z0)  dx0dz0 

XQ  ZO 

+  //  Zzx  (x  —  x0,  z  -  z0)Jx(xo,  z0)dx0dz0  (2.1) 

lo  20 

Ex(x,d,z)  =  J  J  Zxz(x  -  x0,  z  -  z0)Jz(x0,  z0)dxodz0 
*0  20 

+  //  E xx^^  2*  *o)  dx o  dzQ  (2.2) 

Xo  20 

In  order  to  develop  the  integral  equation,  we  enforce  the  boundary  condition 
requiring  that  the  components  of  the  tangential  electric  field,  Ez  and  Ex,  be 
zero  on  the  conductor. 

J  J  Zzz(x  -  x0,z  -  z0)Jz(x0,  Zo)  dx o  dz0 
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+ 


+ 


I  j  Zzx  (a;  —  xq ,  z  —  zq)Jx(xq ,  z0)  dx 0  dz0  =  0  for  x,  z  €  S. 

Xo  20 

^  J"  Zxz ( ^  Xo 5  Z  Zo)  J z  ( Xq,  Zo)  dx o  dzQ 

XO  20 

J  J  Zxx(x  -  Xo,  z  -  zo)  Jx(xq,  zo)  dx0  dz0  =  0  for  x,  z  €  S. 


(2.3) 


(2.4) 


These  equations  constitute  a  set  of  coupled  integral  equations  for  the  unknown 
current  on  the  conductor.  These  two  equations  together  are  equivalent  to  the 
electric  field  integral  equation  (EFIE).  They  are  Fredholm  equations  of  the 
first  type  [24];  that  is,  the  unknown  functions,  J*(xo,  z0)  and  Jx(x0,  z0)  appear 
only  in  the  integrands.  At  this  point  we  can  transform  the  coupled  integral 
equations  into  coupled  algebraic  equations  in  the  spectral  domain  as  in  [26]  by 
making  use  of  the  convolution  theorem  [24]. 


Z„(a,  0)  Jz(a,  0)  +  Z2X(a,  0)JX( a,  0)  =  Ez{a,  0)  (2.5) 

ZX2(a,  0)Jz(a,  0)  +  Zxx{a ,  0)  Jr(a,  0)  =  Ex{a,  0)  (2.6) 


The  following  convention  is  used  for  the  Fourier  transform. 

+oo  +oo 

E(oc,  0)  =  J  J  F(x,  z)  exp(+j0z)  ex-p(+jax)dxdz  (2.7) 

— oo  — oo 

+oo  +oo 

E(x,z)  =  J  J  F(a,0)exp(-j0z)  exp(-jax)  da  d0  (2.8) 

'  *  -oo  -00 

However,  although  this  is  the  usual  method,  it  does  not  make  any  difference  at 
which  point  in  the  analysis  we  transform  to  the  spectral  domain.  Notice  at  this 
point  the  coupled  algebraic  equations  contain  the  unknown  functions  Jz(a,0), 
Jx(a,0),  Ez(a,  0),  and  Ex(a,0).  So  it  appears  we  have  introduced  two  more 
unknown  functions  by  transforming  into  the  spectral  domain.  However,  it  will 
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be  seen  that  the  functions  Ez(a,(3 ),  and  Ex(a,/3)  will  be  eliminated  in  the 
solution  process.  Whether  we  choose  to  solve  the  coupled  integral  equations  in 
the  space  domain  or  the  coupled  algebraic  equations  in  the  spectral  domain, 
we  must  solve  for  unknown  functions.  To  do  this  we  make  use  of  the  method 
of  moments  to  solve  the  coupled  integral  equations  (in  the  space  domain)  or 
coupled  algebraic  equations  (in  the  spectral  domain).  Since  it  is  easier  to 
visualize  the  problem  in  the  space  domain,  we  will  proceed  from  the  coupled 
integral  equations. 

2.2  The  Application  of  the  Method  of  Moments 

As  pointed  out  in  the  previous  section,  the  coupled  space-domain  in¬ 
tegral  equations  and  the  coupled  spectral  domain  algebraic  equations  contain 
unknown  functions.  Therefore,  the  first  step  in  solving  the  equations  is  to  ex¬ 
pand  the  unknown  current  in  an  appropriate  set  of  known  expansion  functions. 

M 

Jz(xq,zo)  =  X]amJ'*po(xo,zo),  and 

1 

M 

•/.(*>,*<>)  = 

1 

where  am  and  bm  are  unknown  constants.  The  required  expansion  functions 
are  dependant  on  the  particular  problem.  They  will  be  discussed  in  detail 
in  Chapter  3.  In  the  deterministic  formulation,  in  addition  to  the  unknown 
current  we  also  have  a  known  source  current  the  x-directed  and  x-directed 
components  of  which  are  given  by 

j;ource(x0,2o)  and 
J*OUTce(xo,  z0)  respectively. 
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Upon  substitution  into  the  EFIE  we  obtain 


J  j  Zzz[x  -x0,z-zQ) 


r  M 


IQ  zo 


E^JlT^o,z0)  +  Jr™(*  o,*o) 


i  1 


dx  o  dzQ 


"b  J  J  ^  ^o) 


XO  zo 


M 


£  bm  J'XT(*0,Zo)  +  JrTCe(xo,Z0) 


L  1 


x,  z  €  5,  and 


dx0  dzo  =  0 
(2.9) 


J  J  Zxz(x  -  x0,  z  -  Zo) 


XQ  Zo 


M 


E  *«)  +  ZO) 


L  1 


c/xq  cLzq 


+  J  J  Zxx{x  -Xo,Z  -  Zo) 


r  m 


Xq 


£  6mCa(^0,  *o)  +  Jrrce(^0,  2o) 


L  1 


x,  z  €  5. 


dx  o  dzo  —  0 
(2.10) 


So  we  have  a  set  of  two  coupled  integral  equations  with  2M  unknown  constants. 
In  order  to  generate  a  linear  system  of  equations,  we  test  the  set  of  equations 
using  an  appropriate  set  of  testing  functions. 


jir  (x,z) . . .  j%\x,z). 
JT{x,z)...J%\x,z). 


As  with  the  expansion  functions,  the  appropriate  testing  functions  depend  on 
the  particular  problem.  These  will  be  discussed  in  detail  in  Chapter  3.  The 
testing  functions  are  used  to  generate  2 N  equations  as  follows. 
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JJjtmJJ 

x  z  ro  zq 

■  |z„(x  -  x„,  z  -  z„)  (E  ■.«*(•  0, zo)  +  Jr“(x«,  z„)) 

+  Z„(x  -  x0, 2  -  Zo)  |^E  bmJ'%a(x 0,  z0)  +  JTr:’txa.  z0)j 


<ixo  dzQ  dx  dz  =  0 


X  Z  To  Zq 

■  [z„(x  -  x„,  z  -  z„)  (E  am  J""“(x0,  Zo)  +  zo)  ) 

+  Z„(x  -  x0,  z  -  z0)  ^E  Zo)  +  J»our«(x0t  zo)) 


<fx0  dzo  dx  dz  =  0 


X  Z  Xq  ZO 

■  ^Z„(x  -  x„,  z  -  zo)  (E  o.^r(»o.  Zo)  +  •/f*r"(x„,  zo)) 

+  Z«(x  -  Xo,  Z  -  Zo)  |E  bmJ'x m“(zO>Zo)  +  7*"ur<'"(x05Zo)) 


6?xq  <i?o  dx  dz  =  0 


X  Z  XQ  Zo 

■  fz„(x  -  Xo,  z  -  z„)  (y.  0,  Zo)  +  rrr“(* 0.  Zo)) 

+  Zrr(x  -  x0,  z  -  Zo)  IE  bmJl%a(x 0,  Zo)  +  7’°"rc*(x o,  Zo)) 


dx0  <izo  dx  dz  =  0 


(2.11) 
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In  the  space  domain,  the  right-hand  side  of  the  integral  equation  is 
zero  because  the  electric  field  is  zero  on  the  conductor.  However,  the  corre¬ 
sponding  spectral  domain  algebraic  equation  is  not  zero  because  the  Fourier 
transform  covers  the  entire  y  =  d  plane.  In  the  eigenvalue  approach  [26],  Par- 
seval’s  theorem  is  used  to  show  that  the  right-hand  side  of  the  tested  equations 
is  zero.  The  same  situation  exists  for  the  excitation  case.  Note  here  that  the 
right-hand  side  of  the  tested  space-domain  equations  is  zero.  If  we  had  teoted 
the  equations  in  the  spectral  domain,  the  right-hand  side  would  still  be  zero 
by  Parseval’s  theorem.  The  source  term  is  embedded  in  the  left-hand  side  of 
the  equation.  In  some  situations  it  will  be  desirable  to  use  a  different  number 
of  testing  functions  than  expansion  functions  to  generate  an  overdetermined 
system  of  equations.  In  the  case  of  an  excitation  problem,  in  general,  it  is  no1- 
possible  to  make  use  of  Galerkin’s  method1. 

The  tested  equations  result  in  a  linear  system 


'  [KZZ][KZX)  ' 

to] 

[  IS.) 

.  [KXZ}[I<XX] 

[  l*]  J 

L  [£*] 

where  [a]  and  [6]  are  vectors  of  length  M  containing  the  unknown  constants 
a\ . . .  a.M  and  fej  . . .  b\f  in  the  expansions  for  Jx  and  Jz ,  [ Kzz\  ■ . .  [ Kxx\  are  M  x  .V 
matrices  given  in  the  space  domain  by 

Actually,  while  it  is  not  possible  to  apply  Galerkin’s  method  to  a  deterministic  problem 
using  a  travelling  wave  source  formulation,  it  is  possible  to  do  so  .  hen  a  localized  source  is 
used.  For  example  a  planar  dipole  driven  by  a  current  source  at  the  center  could  be  analyzed 
using  Galerkin’s  method  because  the  current  exists  only  over  a  finite  area. 
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4-oo  +oo  4-oo  +oo 


A',7  =  J  J  J  J  Jt'st{x,z)Zzz(x-  Xo,Z-  Zo) 


—  OO  —  OO  — OO  — OO 


•  JzjPa(x0,  zo)  dx  dz  dx o  dz0 , 


+oo  +oo  4-oo  +oo 


A',7  =  J  J  J  J  Jtzeiat(x,z)Zzx(x-x0,z-  Zo) 


—  OO  —co  — OO  — OO 


•J'*pa(x o,  z0)  dx  dz  dx o  dz0, 


4-oo  4-00  4-oo  4-00 


=  /  //  J  J«‘(x,z)Z„(x-x0,z-z0) 


— oo  — oo  — oo  — oo 


■Jl*pa(x o,  zo)  dx  dz  dx o  dzo,  and 


4-oo  4-oo  +oo  4-oo 


A'7  =  J  J  J  j  2)Zxx(x  -X0,Z-  Zo) 


—  oo  — oo  — oo  — oo 


■J‘*pa(xo,  zo)  dx  dz  dx o  dz0 


and  [Sz]  and  [Sr]  are  vectors  of  length  M  given  by 


+oo  +oo  +oo  +oo 


S’  =  -  j  j  j  J  Jf‘(x,z)Z2z(x-x0,z-z0) 


-oo  — oo  — oo  -oo 


■j;ourc'(Xo,zo)dxdz  dxodzo 


+oo  4-oo  4-oo  4-oo 


~  J  J  J  J  dleiat(x^z)Zzx{x  ~  X0,Z  ~  Z0) 


— oo  — oo  — oo  — oo 


■J’OUTCe(x0,  z0)  dx  dz  dx0  dz0 ,  and 


4-oo  4-oo  4-oo  4-oo 


Sf  =  -  j  f  I  J  JTMzxt{x-x0,z-zo) 


-OO  — _ _  — OO  —  OO 


j  source  ^  fa  fa  faQ  faQ 


4-co  4-oo  4-oo  4-oo 


-  J  J  J  J  Jt'*\x,z)Zxx{x  -  x0,  z  -  Zo) 


—  OO  — OO  — OO  — OO 


•  J*OUTCe(xQ ,  Zo)  dx  dz  dxo  dz0. 


(2.12) 

(2.13) 

(2.14) 

(2.15) 

(2.16) 

(2.17) 


However,  since  the  Green’s  function  is  only  available  in  closed  form  in  the 
spectral  domain,  these  integrals  are  six-dimensional  and  therefore  extremely 
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expensive  to  evaluate.  Some  work  has  been  performed  to  reduce  this  complexity 
[16].  In  order  to  transform  the  integrals  into  the  spectral  domain,  we  express  the 
space  domain  Green’s  function  as  the  inverse  Fourier  transform  of  its  spectral 
domain  counterpart.  For  example 

4-oo  +oo 

Zzz(x  -  x0,z  -  z0)  =  J  J  Zzz(a,/3)exp(-jt3{x  -  x0)) 

— oo  — oo 

•  exp(— ja(z  —  z0))dad/3  (2.18) 

We  then  substitute  this  expression  into  the  integral  to  obtain 

+oo  -foo  +oo  +oo 

K$  =  /  /  j  I 

— oo  — oo  — oo  — oo 
-foo  +oo 

•  J  J  Zzz(a,  /3)  exp(—j/3(x  -  x0))  exp(-ja(z  -  z0))dad(3 

— oo  — oo 

■Je2Jpa(x0,z0)dxdzdxodz0.  (2.19) 

Rearranging  the  order  of  integration,  we  obtain 

+oo  +oo  r  +oo  4-oo 

K**  —  f  f  f  j  z)  exP(~j @x)  exv(-jQZ)dx  dz 

— oo  — oo  L-oo  — oo 

•  Z„(a,0) 

4-oo  -foo 

//  Jz*pa(x0,  zq)  exp(j0xo)  exp(jaz0)dx0  dz0\  dadfd.  (2.20) 

-oo  — oo 

We  recognize  that 

4-oo  4-oo 

jt*atm(ai,  0)  =  J  J  J‘“‘(x,  z)  exp(— ;/?x)  exp(—jaz)dx  dz,  and  (2.21) 

-OO  -oo 
-foo  -foo 

J"pa(a,/ 3)  =  J  J  Jez*pa(xo,z0)exp(j@x0)exp(jazQ)dxodzo.  (2.22) 

-oo  -oo 
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Therefore,  in  the 

spectral  domain,  [ Kzz\ . . .  [Kxx]  are  given  by2 

Kzz(i,j) 

+oo  +oo 

(2.23) 

Kzx(i,j) 

— oo  —  oo 

+oo  +oo 

(2.24) 

Kxz(i,j) 

“OO  “OO 

+oo  +oo 

(2.25) 

K  xx(ij) 

“OO  “OO 

+oo  +  0O 

“OO  — oo 

(2.26) 

and  [Sz]  and  [Sz] 

are  vectors  of  length  M  given  by 

si  = 

+oo  +oo 

1  J  H)Z„(a,l3)jr™(a,l3)  da  J0 

+ 

“OO  “OO 

+oo  +oo 

J  J  jir-la, /3)Z«(a,0)jr'a(a,l3)  da 

(2.27) 

Sf  = 

— oo  — oo 

+oo  +oo 

J  J  Jg*(a, 0)Z„(a,l))jrr‘’(a,l3)dadl3 

+ 

“OO  “OO 

+00  +oo 

J  J  JlV(a,l3)ZII(a,l3)jrr“(a,0)dad0. 

“OO  — OO 

(2.28) 

The  linear  system  is  now  solved  for  a\ . . .  a\f  and  bx . .  .b\f  using  Cans- 

sian  elimination  if  the  system  is  square.  If  the  system  is  overdetermined  we 


generate 


'  [K„][K„\  ' 

T 

'  [*„]  \i<„\ ' 

[«]  1 

■  [*..]  l K„]  ' 

T 

f  IS.] ' 

.  [K„\  [Kxx] 

.  {K„){K.A  . 

[b}\ 

.  [*„][*„] 

[  IS*)  J 

2Notice  that,  in  the  formulation  given  here,  the  convention  for  the  Fourier  transform  is  the 
opposite  of  that  used  in  [15].  Therefore,  the  conjugation  in  the  inner  products  is  performed 
on  the  testing  functions  as  opposed  to  the  expansion  functions. 
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where  the  superscript  “T”  denotes  the  transpose  of  a  matrix.  We  can  then 
solve  the  resultant  2M  x  2 M  linear  system  of  equations.  This  is  equivalent  to 
a  least  squares  formulation. 


2.3  The  Green’s  Function 

The  characteristics  of  the  Green’s  function  are  important  because 
much  information  may  be  derived  from  them.  The  spectral  domain  Green’s 
functions  used  in  the  following  analysis  are  derived  using  the  immittance  ap¬ 
proach  [26]  in  Appendices  1  and  2. 

The  spectral  domain  Green’s  function  for  the  grounded  dielectric  slab 
is  given  by 

Zzz  =  - N2zZ'-N2xZh ,  (2.29) 

Zzx  =  -NzNx(Ze  —  Zh),  (2.30) 

Zxz  =  —NzNx(Ze  —  Zh)  and,  (2.31) 

Zxx  =  —NlZ*  -  N]Zh,  (2.32) 


where 


ne  =  Ytmu 


(2.33) 

(2.34) 

(2.35) 

(2.36) 

(2.37) 
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Y2e 

= 

Yr  M2  coth(72d), 

(2.38) 

ff 

= 

Ytei, 

(2.39) 

K* 

A. 

— 

YtE2  rot.h(72d), 

(2.40) 

Ytm\ 

= 

jve  o 

> 

7i 

(2.41) 

Ytm 2 

= 

jutRdo 

•) 

72 

(2.42) 

Ytei 

= 

7i 

jojp' 

(2.43) 

YtE2 

= 

72 

jup' 

(2.44) 

7i 

= 

\Ja2  +  /32  —  k2  and, 

(2.45) 

72 

= 

Ja2  +  /32  -  eRk2. 

(2.46) 

It  can  be  seen  that  the  components  of  the  Green’s  function  have  poles  corre¬ 
sponding  to  the  poles  of  Ze  and  Zh.  The  poles  of  Zc  correspond  to  the  TM 
surface  waves,  and  the  poles  of  Zh  correspond  to  the  TE  surface  waves.  The 
poles  occur  when 


e/Oi  +  72  tanh(72d) 
7i  +  72  coth(72d) 


0,  and 

0. 


(2.47) 

(2.48) 


In  order  to  facilitate  the  study  of  the  surface  wave  poles,  we  perform  a  change 
of  variables  defined  by 


p  =  \jo. 2  +  /?2,  and  (2.49) 

6  =  tan-1  .  (2.50) 

\  r  / 


We  note  that  Ze  and  Zh  are  functions  of  p  only.  In  this  polar  coordinate 
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Region 

Range 

7i 

72 

1 

p  <  k 

imaginary 

imaginary 

2 

k  <  p  <  y/e^k 

real 

imaginary 

3 

\/£Rk  <  p 

real 

real 

Tablf*  2.1:  Regions  of  the  Radial  Wavenumber,  p 


system,  the  integrals  take  the  following  form. 

+oo  2jr 

=  f  j (2.51) 
0  0 


Now  we  define  three  regions  as  in  Table  2.3.  All  of  the  surface  wave  poles  occur 
in  region  2.  This  is  occurs  because,  in  order  for  a  surface  wave  mode  to  exist, 
the  field  outside  the  slab  must  be  evanescent  in  the  y-direction  while  the  field 
inside  must  be  a  standing  wave  in  the  y-direction.  The  values  of  p  at  which  the 
surface  wave  poles  occur  are  the  values  of  the  propagation  constants  for  the 
surface  wave  modes.  In  order  to  evaluate  the  integrals,  we  must  circumvent 
the  poles.  The  contributions  of  the  poles  are  included  via  the  calculation  of  the 
residues  at  the  poles.  We  can  show  that  the  poles  are  of  the  first  order.  Using 
the  integration  path  given  in  Figure  2.2  and  the  Cauchy’s  residue  theorem,  we 
can  show  that  when  only  the  T Mo  mode  can  propagate,  Kzz(i,j)  ■  ■  ■  Kx x(i,j)  are 
given  by 

PTM0~&  2tt 

=  J  J j?“lp,<>)Z"(p,e)J;r(p,i>)pdpde 
0  0 

-  iitRes  j  (J?'-(p,<))Z„(P,<>)j;r(P’l»)  \^rmd9 
0 
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+oo  2jt 

+  /  j  <2-52) 

PTMO+&  0 


The  singularity  occurs  only  in  the  Green’s  function,  so  we  write 

Ptmo-6  2ir 

K„0.i)  =  /  j  3T-(e,<>)Z,.(p,<>)jl7'(p,9)dadl3 

0  0 
2ir 

-  it  j  J,'?“(p,e)Z„TM(p,6)j“’°(p,l))dadl3 

0 

oo  2ir 

+  /  J j?‘-(p,e)z,2(p,e)j;*’”(pj)dadi3,  (2.53) 

PTMO+5  0 
PTMO-S  2ir 

=  /  J  jlf(p,e)za(p,e)j"’’‘(p,o)dad0 

0  0 
2ir 

-  i*  J  JT'MZ'-TuipJVlT'Mdadt} 

0 

00  2ir 

+  J  j  jT’MZ„(p,9)j^{p,e)dcd3,  (2.54) 

ptmo+6  0 
PTMO-5  2tt 

=  /  j  jT'MZ*'MjT^e)dadl3 

0  0 
2* 

-  in  j  J'sf'ip,  0)ZtzTM{p,  9)Jez*pa{p,  9)  doc  dp 

0 

00  2tt 

+  /  j  J,"“(p,»)ZI,(p,e)j"’n(pJ)dQd0,!i  nd  (2.55) 

PTMO+&  0 
PTMO-t  2  if 

A-,,!,.,)  =  /  / JT‘(P’t>)z**(p>0)J"’°{p,e)dadt) 

0  0 
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2ir 

-  IX / J‘,f(P,O)Z,ITM(pJ)J;7'‘(p,e)dad0 
0 

oo  2tt 

+  /  //“"■(/-,  S)Z„(pJ)j"^(p,O)dad0  (2.56) 

PTMO+5  o 

where 


7e 

^ TM 

dY‘ 

dp 

dJl 

dp 

For  T E  modes  the 


ZzzTM 

=  ~N2zZeTM, 

(2.57) 

ZzxTM 

=  - NzNxZ\m , 

(2.58) 

ZxzTM 

=  -N2NxZ'tm, 

(2.59) 

ZxxTM 

(2.60) 

1 

(2.61) 

1 

i 

=  (^)n*w.^d 

(2.62) 

=  (^r)  (*?(/>)  +  V-TM2rf72(cothJ(7,c()  -  1)) 

(2.63) 

residue  is  given  by 

ZzzTE  = 

—  /V2  7^ 

iViZT£t 

(2.64) 

ZZxTE  = 

-\-NzNxZjE, 

(2.65) 

^izTE  = 

+NiNxZ2e  and, 

(2.66) 

ZxxTE  ~ 

-NlZiE 

(2.67) 

where 


1 

dYp  3Y£_  ’ 
9p  '  3p 


(2.68) 
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=  ^  and  (2.69) 

=  U)(Y^)  +  YTE2d^M-l)).  (2.70) 

It  should  also  be  noted  that  the  Green’s  function  has  a  branch  point  at  p  =  k0. 
However,  the  Green’s  function  is  bounded  at  the  branch  point  and  therefore, 
no  numerical  difficulty  is  encountered  except  the  correct  choice  of  Riemann 
sheets.  This  amounts  to  choosing  71  to  be  positive  real  for  p  >  kQ  and  positive 
imaginary  for  p  <  ko.  There  is  no  branch  point  at  p  =  ^/eriko-  This  is  because 
72  always  appears  with  coth(72<i)  in  the  Green’s  function  thus  cancelling  the 
sign  of  72.  Therefore,  the  sign  of  72  need  not  be  specified.  The  pole  in  the 
Green’s  function  corresponding  to  a  particular  surface-wave  mode  occurs  at 
the  branch  point  at  the  cutoff  frequency  of  the  mode  and  then  migrates  toward 
tnko  with  increasing  frequency.  Therefore,  the  surface  wave  poles  occur  in  the 
order 


dYf 

dp 

dY2h 

dp 


k0  <  .  .  .  PTE3  <  PTM2  <  PTEl  <  PTMO  <  y/^R^o 


(2.71) 


2.4  Eigenvalue  Formulation 

The  eigenvalue  formulation  is  used  here  to  analyze  a  uniform  line.  We 
assume  that  the  current  on  the  line  is  of  the  form 


Jz(x,z)  =  Jz(x)  exp(—j  f3zz) 
Jx{x,z )  =  Jx(x)exp(-j0zz) 


(2.72) 

(2.73) 
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lm(p) 


Figure  2.2:  Integration  Path 


We  make  use  of  a  Galerkin  method,  that  is  we  make  the  testing  functions 
identical  to  the  expansion  functions.  Then  we  generate  the  determinant  of  the 
coefficient  matrix.  This  determinant  is  a  function  of  (3.  We  determine  (3  =  (3Z 
such  that  the  determinant  is  zero.  This  is  equivalent  to  solving  the  nonlinear 
equation  for  (3 


[Kzz]  [Kzx] 
[K„]  [K„] 


=  0. 


This  equation  is  solved  using  an  interval  halving  technique.  Once  f3z  is  known, 
we  determine  the  transverse  dependance  of  the  current  from  the  eigenvectors. 
The  calculation  of  the  characteristic  impedance  is  not  unique  because  the  fun¬ 
damental  mode  is  not  perfectly  TEM.  Two  definitions  are  possible. 


Z0 


(2.74) 


where  Pz  is  the  z-directed  power  flow,  Iz  is  the  z-directed  current  flow,  and 
V  is  the  potential  between  the  strip  and  the  ground  plane.  We  use  the  first 
definition  because  V  is  not  unique.  The  z-directed  current  on  the  strip  is  given 
by 

w/2 

Jz=  Jz(x)dx.  (2.76) 

-m/2 


The  z-directed  power  flow  is  given  by 


Pz  =  Re 


OO  OO 

J  J  (E(x,y)  x  H*(x,t/))  •  zdxdy 


.0  —  oo 


(2.77) 


The  x-portion  of  the  integration  may  be  transformed  into  the  spectral  domain 
through  the  use  of  Parseval’s  theorem. 


Pz  =  —Re 


/M  a,y)xH-(a,y))-zdady 

.0  —  oo 
JO  oo 

=  2^Re/  /  [(£*1  {Q,y)Hh{Q,y)- 

d  —  oo 

Eyi(°‘,y)Hxi{<*,y))da  dy ] 
d  oo 

+  2tt  Re/  /  [(Ex2{a,y)Hl2(a,y)- 

0  —oo 

Ey2(a,y)H*2(a,y))da  dy] 


(2.78) 


(2.79) 


where 
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Hxi{a,y) 

= 

exp(-7iy), 

Hyi  (a,y) 

= 

^Ji(a)exp(-7iy), 

£*2(0,2/) 

= 

£'°2(Q)sinh(72y), 

Ey2(ot,y) 

= 

£j2(a)c°sh(72y), 

HX2  (ct,y) 

= 

#°2(o)  cosh(72y), 

Hy2  (o,y) 

= 

#J2(a)  sinh(72y), 

&(«) 

= 

-r(?-)Be  +  9£A, 

y\ 

£» 

= 

«?,(<*)  =  «B.  +  r(?-)B», 

2l 

K(°)  = 

i'”2(a)  =  ryA,  +  «/4», 

£»  =  f-4., 

S/2 

=  -?Ae  +  r^, 
^2 

BS,(a)  =  (f)A», 

22 

a 

r  =  -r, 

P2 

P 

q  P2' 

V\  =  jweo, 

j?2  =  jwe0e«, 

2i  =  jw/i, 


2  2 


jW/i, 
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A,  =  ~l(a)p, 

(-'i  Vi 

Ah  = 

B,  =  D^j,(a)p, 

yi 

Bh  =  dAJMp, 

7i 

Ci  =  T^sinh(72tf)  +  cosh(72rf), 

V2  y  i 

C2  =  sinh(72rf)  +  — cosh(72</), 

7i  22 

n  exp(7i^)sinh(72cf) 

Dl  =  - ci - 1 

„  exp(7i£/)sinh(72<f) 

.  =  ~ - 
ju{a)  =  NxJz(a)  -  NzJx(a),  and 

J.(o)  =  N,J,(a)  +  Ngjx{a). 

The  y-integration  is  simple  and  may  be  done  in  closed  form.  Therefore,  we 
may  write 

{Si{a)  +  S2{a))da  (2.80) 


where 

S,  =  (jM  (£»  (o)WjT(a)  -  E*  (a)/ft(a))  exp(-2 7i<0,  and  (2.81) 
S2  =  [^](£’°2(Q)W°2'(a)sgn(7j)(— sinh(7jii)cosh(72d)  -  d] 

*  72 

-  E%(a)H°x;(a)[-smh(^d)cosh(72d)  +  d}].  (2.82) 

72 


where  sgn  is  1  when  72  is  real  and  -1  when  72  is  imaginary.  The  a  integration 
is  done  numerically. 
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2.5  Source  Formulation 

The  source  formulation  is  the  crucial  part  of  this  method.  The  source 
formulation  affects  the  accuracy  and  efficiency  of  the  method  and  the  ease 
of  determining  network  parameters.  Several  source  formulations  have  been 
proposed.  The  one  used  here  appears  to  have  first  been  used  by  [17].  It  was 
later,  perhaps  independently,  proposed  by  [15].  This  method  involves  the  use 
of  a  traveling-wave  formulation.  First,  a  spectral  domain  eigenvalue  method 
is  used  to  determine  the  propagation  constant  and  current  distribut  on  for  a 
uniform  line.  Then  this  information  is  used  to  derive  the  source  function,  wh:ch 
consists  of  an  incident  current  wave. 

jrTC'(x,z)  =  £>,  ./,,(*)  exp  (-j0z)  (2.83) 

i=i 

•/“"•(l,*)  =  X>.4,(*)exp(-j7Jz)  (2.84) 

1=1 

The  functions  Jzm  and  Jxn  are  the  Maxwellian  basis  functions  which  will  be 
discussed  in  the  next  chapter.  They  include  the  proper  edge  condition  for 
the  longitudinal  and  transverse  current  at  the  edges  of  the  microstrip.  The 
constants  ax . .  .  ai  ,  bx  . . .  6^,  and  /?,  the  propagation  constant,  are  determined 
in  the  spectral  domain  eigenvalue  analysis. 

Although  in  scattering  formalism  the  source  is  postulated  as  an  in¬ 
coming  wave  from  2  =  —00,  the  source  function  is  truncated  at  some  point  in 
oxder  to  facilitate  the  numerical  integration  of  the  inner  products.3  We  will  now 

3It  is  not  necessary  that  the  source  functions  be  truncated  in  order  for  the  innner  products 
to  exist.  This  will  be  discussed  in  Chapter  3. 
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show  why  the  integrals  may  be  truncated  while  maintaining  the  correctness  of 
the  formulation. 

Consider  the  electric  field  generated  by  a  traveling  wave  on  an  infinite 
line.  For  simplicity,  consider  only  longitudinal  current  on  the  line. 

Ez  -  J  J  Zzz(x  -  x0,z  -  zo)^2aiJzi(xo)exp(-j/3zzo)dxodso  (2.85) 


In  the  spectral  domain  this  becomes 


Ez{a,/3)  =  Zzz(a,l3)Y2aiJzm(ct)6(f3z). 


(2.86) 


Therefore,  after  testing  we  have, 


+oo  +oo 


J  J  j‘r,(a,0)Z„(aJ)jf’'-(a,l3)dad0 

-oo  -oo 
+oo 


(2.87) 


by  the  sifting  theorem.  In  the  eigenvalue  formulation  we  determined  (3  such 


“TOO 

J  j‘;-,(a)Z„(a,0,)jt‘”"‘(a)dc.  =  0. 


(2.S8) 


So  we  see  that  a  traveling  wave  on  an  infinite  line  automatically  satisfies  the 
integral  equation  for  the  electric  field.  Likewise,  a  standing  wave  on  an  infinite 
line  satisfies  the  same  conditions.  It  seems  reasonable  to  assume  that  a  traveling 
wave  on  a  finite  line  satisfies  this  condition  far  away  from  any  discontinuities. 
Therefore,  if  we  consider  an  incident  wave  on  a  semi-infinite  microstrip  line 
from  — oo,  the  electric  field  is  given  by 
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Ez(x,d,z) 


Ex{x,d,z) 


+w  0 


—  J  J  Zzz(x  Xq ,2 

—w  — oo 
+u;  0 

"H  J  J  Zzx(x  x0,  z 

—  til  — oo 
+«/  ~ZL 

=  j  j  ^zz(%  ^0,  z 

—  u/  — oo 
+ti>  ~ZL 

+  J  J  Zzx{x-x0,z 

—  til  —CO 

w  0 

■f"  j  j  ^ zz ^  ®0,  ^ 

-tu  -zl 
w  0 

+  J  J  Zzx(x-x0,z 

-tii  -Zi 
+  U;  0 

~  J  J  Zxz(x  —  xo,  z 

— w  — OO 
+tx/  0 

+  II  ^ xi (*^  *^0i  ^ 

— ty  — oo 
+w  ~zL 

~  J  J  Zxz(x  xQ,z 

—w  — oo 

+  W  -  ZL 

+  J  J  Zxx(x-X0,z 

—w  —  oo 
+u/  0 

+  J  J  Zxz(x-x0,z 

-w  -  Zl 
+w  0 

-h  J  J  Zxx(x  x0,  z 


-  zo)Jz(xo,  20)  dx o  dz0 

-  z0)jx(x0,  z0 )  dx o  dz0 
20)J z(^o,  2o)  dx o  dzQ 

-  Zq)Jx(xq ,  20)  ^0 

-  20)  J2(x0,  20)  <fc0  d20 

-  zo)./x(zo,  zo)  c?xo  dzo,  and  (2.89) 

-  Zo)Jz(x0,  z0)  dx o  dz0 

-  Zo)Jx{x0,  z0)  dx o  dz0 
zo)Jz(x0,z0)dx0  dz0 

~  zo)Jx(x0,  zq)  dx o  dz0 

-  zq)Jz[xq,  Zo)  dx o  dzo 

-  z0)Jx(x0,  z0)  dx0dz0.  (2.90) 


-w  -  zl 


If  2  >>  —  zl,  then  the  integrals  extending  from  — oo  to  —  zl  approach  zero. 
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Therefore,  it  is  equivalent  to  define  the  Fourier  transform  of  the  travelling 
wave  current  as 

+oo  0 

a  J(x,z)  exp(jax )  exp (jf3z)dx  dz.  (2.91) 

-CO  -ZL 

This  analysis  may  be  extended  to  the  two-dimensional  case  simply  by  noting 
that  the  uniform  line  current  is  the  homogeneous  solution.  Therefore,  it  may 
be  added  to  the  unknown  current  without  affecting  the  problem. 

The  reflected  wave  and,  in  the  case  of  a  multiport,  the  transmitted 
wave(s)  may  be  treated  in  a  similar  manner.  The  point  at  which  the  functions 
are  truncated  affects  the  rapidity  of  the  convergence  of  the  integrals.  This  will 
be  discussed  in  the  next  chapter. 

The  equation  is  tested  in  order  to  enforce  the  boundary  condition, 
Eton  =  0  on  the  strip.  The  test  region  must  extend  far  enough  from  the 
discontinuity  that  scattered  field  is  taken  into  account.  However,  it  must  not 
extend  too  close  to  the  point  at  which  the  uniform  functions  are  truncated. 


Chapter  3 


DETAILS  OF  THE  NUMERICAL  COMPUTATION 


This  chapter  deals  with  the  details  of  the  computation.  Because  full- 
wave  methods  such  as  the  one  presented  here  are  computationally  intensive,  it 
is  necessary  to  take  advantage  of  all  available  numerical  labor-saving  devices. 
Therefore,  we  will  discuss  the  following  aspects  of  the  computation. 

•  The  choice  of  expansion  and  testing  functions; 

•  The  use  of  symmetry  properties  in  the  reduction  of  the  numerical  inte¬ 
gration; 

•  The  use  of  symmetry  properties  in  filling  of  the  impedance  matrix; 

•  The  conditioning  of  the  integrands  through  the  extraction  of  the  singu¬ 
larities  associated  with  the  surface-wave  poles; 

•  The  truncation  of  the  traveling  wave  functions  for  optimum  convergence; 

•  The  acceleration  of  convergence  of  the  integrals  through  the  extraction 
of  asymptotic  forms; 

•  The  details  of  the  computational  algorithm. 
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3.1  Expansion  and  Testing  Functions 

The  number  of  basis  functions  required  may  be  minimized  if  basis 
functions  represent  the  current  on  the  conductors  as  closely  as  possible.  The 
computational  effort  may  thus  be  reduced  significantly.  The  basis  functions 
should  exhibit  the  proper  edge  singularities,  boundary  conditions,  and  sym¬ 
metry.  Since  the  computations  are  to  be  performed  in  the  spectral  domain, 
it  is  desirable  to  be  able  to  obtain  the  Fourier  transforms  in  closed  form. 
Both  entire  domain  and  subsectional  basis  functions  are  used  in  this  study. 


Two-dimensional  basis  functions  are  expressed  as  separable  products 

dimensional  function  of  x  and  a  one  dimensional  function  of  z. 

of  a  one 

Jtxpa(x,z ) 

(3.1) 

Jlest(x,z) 

(3.2) 

Jexxpa(x,z) 

=  r^(x)rx^(z) 

(3.3) 

J?*(*,z) 

=  JS3(*)^33(») 

(3.4) 

Therefore,  the  Fourier  transforms 

are  expressed  as  separable  products 

of  cv  and 

j:xpa(c*,P) 

=  nW 

(3.5) 

jj-W) 

=  JSZWJSi )W 

(3.6) 

jlxpa{a,  (3) 

- 

(3.7) 

(3.8) 

Figure  3.1:  The  Maxwellian  Functions 


3.1.1  Analysis  of  Uniform  Microstrip  Line 


In  the  one-dimensional  analysis  of  uniform  microstrip  line,  entire  do¬ 
main  basis  functions  and  Galerkin’s  method  were  used.  The  Maxwellian  basis 
functions  were  chosen  to  represent  the  transverse  dependance  of  the  current 
on  the  uniform  portions  of  the  conductors.  These  functions  exhibit  the  proper 
symmetry  and  accurately  include  the  effects  of  edge  conditions. 


(3-9) 


•J arm  (  &  ) 


(3.10) 


The  one-dimensional  Maxwellian  basis  functions  are  shown  in  Figure  3.1.  This 
choice  has  been  shown  [26]  to  give  very  accurate  results  with  a  minimal  number 
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of  functions  due  to  the  variational  formulation. 


3.1.2  Analysis  of  Discontinuities 

When  analyzing  discontinuities,  two-dimensional  basis  functions  in 
which  entire  domain  functions  were  used  to  describe  the  transverse  depen- 
dance  of  the  current  and  subsectional  basis  functions  were  used  to  describe  the 
longitudinal  dependance.  The  Maxwellian  functions  described  above  were  used 
for  the  transverse  dependance.  Several  types  of  subdomain  functions  were  used 
for  the  longitudinal  dependance.  It  was  found  that  piecewise  linear  (overlap¬ 
ping  rooftop)  functions  allowed  approximate  satisfaction  of  all  types  of  edge 
conditions  encountered.  These  functions  are  shown  in  Figure  3.2. 

*«.<*>  =  p-1;-*-1)  (3.H) 


Jzmn(x,z)  = 


=  (-  2nl 


■  h  '  \A  -  (¥ )2 

sinl2^) 


(3.12) 

(3.13) 


■  h  >  \/l-(f)2 

For  example,  the  two-dimensional  basis  functions,  Jz i„(x,  z )  and  Jx in(.r,  z),  are 
shown  in  Figure  3.3  and  Figure  3.4. 

Travelling  and  standing  waves  were  modelled  by  combining  two  two- 
dimensional  basis  functions  in  which  the  Maxwellian  functions  were  used  for  the 
transverse  dependance  and  sinusoidal  functions  for  longitudinal  dependance. 


Jz»in{x,z)  =  sin  (/3zz) 


cos^-ni?) 

i/i  -  (V)a 


(3.14) 


37 


Figure  3.4:  The  2-dimensional  Basis  Function,  Jx\n(x,z) 

(3.15) 

For  example,  a  positive  travelling  wave  can  be  made  up  of  two  phase-shifted 
sine/Maxwellian  functions. 

exp (j/3zz)  =  cos (0zz)  -  j  sin (0zz) 

7T 

=  sin(/3z2  sin(/?2z) 

The  truncation  of  the  sine  function  is  detailed  in  section  3.5. 

3.2  Symmetry  Properties  of  the  Green’s  Function 

The  spectral  domain  Green’s  function  has  the  symmetry  properties 
given  in  Table  3.2.  Because  all  of  the  expansion  and  testing  functions  used  are 


(3.16) 

(3.17) 


—  sin(/322) 
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Component 

Symmetry 

Zz,{a,P) 

even  wrt  a  wrt  ft 

Zzx(a,p) 

odd  wrt  a  wrt  p 

Zxz(a,P) 

odd  wrt  a  wrt  P 

Zxx{ot,  P) 

even  wrt  a  wrt  P 

Table  3.1:  Symmetry  Properties  of  the  Green’s  Function 

purely  real  in  the  space  domain,  the  real  parts  of  their  Fourier  transforms  are 
even  functions  and  the  imaginary  parts  are  odd.  Because  of  these  symmetry 
properties,  we  can  reduce  the  range  of  integration.  The  real  parts  of  the  Fourier 
transforms  of  the  basis  functions  are  always  even,  and  the  imaginary  parts  are 
always  odd.  Thus  the  integration  reduces  to 

+oo  +oo 

*■„( =  J  J 

— oo  — oo 
+oo  +oo 

=  j  j  Z»(a.« 

o  o 

(Re  [iSiW]  «*[•%>(“)] 

+  Imfe(Q)j  Im  [./£»]) 

(Re  [./$&(/?)]  Re 

+  Im  Im  [^J(«])  dccdS,  (3.18) 

+oo  -fco 

R«( =  /  J  J5*(a,0)Z„(a,0)j%~{a,0)dadl} 

— oo  — oo 
■f  oo  -f*oo 

=  j  j  Z-(a,/») 

0  0 

(Re  [.%(<*)]  Im  fe(a)] 

-Imp^lallRe^yc,)]) 


39 


I< 


xz(i,j) 


K 


xx(i,j) 


(Re  [ JS&m]  Im  [jltmW] 

-  Im  [ J$%A0)]  Re  [j^S)iW})  •‘“W, 

+oo  +oo 

j  J  j‘f(a,/3)Z„(Q,/3)j“PV,/9)*<<;/3 

-oo  -oo 


■f  OO  -f  OO 

j  j  Z„(a,ff) 

0  0 

(Re  [./£',;(<*)]  Im  [j£j(«)] 
-Im[j^(a)]Re[J^(a)]) 

(Re  [.%(/»)]  Im  [j^w] 

-  Im[J5&(j9)]lu[jJg£(/J)])  da  iff,  and 

•foo  +oo 

J  J  J*‘at(a,  ff)ZIX(a, ff)J'*pa(a,  ff)dadff 


+oo  +  °° 

/  /  Zxx(a,ff) 

0  0 

(Re  few]  Re  te(«)] 

+  Im  few]  Im  [%(«)]) 

(Re  [Jy&W]  Re  [J^,W] 

+  Im  [j^W]  Im  \j"^(ff)  ])  da  iff, 


(3.19) 


(3.20) 


(3.21) 


3.3  Symmetry  Properties  of  the  Matrix 

A  large  savings  of  computational  effort  can  be  obtained  by  taking 
advantage  of  the  symmetry  properties  of  the  matrix.  In  general,  the  portion  of 
the  matrix  generated  with  Galerkin’s  method  applied  to  the  identical  shifted 
subdomain  expansion  functions  is  Toeplitz  or  anti-Toeplitz.  That  is  the  matrix 


40 


is  of  the  form 

/  7\  T2  T3  T4  \ 

T-  2  Ti  r2  r3 

r_3  t_4  r2 

V  t_4  t_3  r_2  t\ 

This  symmetry  exists  because  the  elements  depend  only  on  the  difference  of 
centers  of  the  testing  and  expansion  functions.  For  example 

+oo  +oo 

KZz(i,i)  =  J  J  ZZz{aiP) 

0  0 

(-%.•(«)%(«) 

Jz(0) o(P)c°s  (0 latest  ~  Zexpa I)  da  dfd  (3.22) 

+  OO  +03 

Kzx(i,j)  —  J  J  Zzx(<X,P) 

0  0 

(*%<(«)*%•(«) 

Jl(t3)o(0)sin  (Potest  -  Zexpa  ))  dad/3  (3.23) 

+O0  +00 

Kxz(i,j)  =  J  J  zxz(a,0) 

0  0 

(P)sin(P(zi',t-  Z'xr*))  dadfi  (3.24) 

+  0O  +OO 

Kxx(i,j)  —  j  j  ^ XX ^ 

0  0 

dlw o(/3)c»*  -  Zerpol)  <Jq  J/)  (3.25) 

Therefore,  for  Kzz  and  Kxx,  T_n  =  Tn  (Toeplitz).  For  Kzx  and  I\xz,  T-n  =  —  Tn 
(Anti-Toeplitz).  In  general,  portions  of  the  matrix  which  are  not  Toeplitz  are 
still  symmetric  about  the  principle  diagonal  due  to  reciprocity.  Furthermore, 
the  zx  and  xz  matrices  are  identical  due  to  reciprocity. 
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3.4  Surface  Wave  Poles 

As  was  pointed  out  in  Chapter  2,  the  spectral  domain  Green’s  function 
has  poles  corresponding  to  TM  and  TE  surface  waves.  These  poles  are  the 
solution  to  the  nonlinear  equations 

cr7i  +  72  tanh(72d)  =  0,  and  (3.26) 

7i  +  72  coth(72d)  =  0.  (3.27) 

These  equations  are  solved  using  an  interval-halving  routine.  The  location  of 
the  surface-wave  poles  must  be  known  accurately.  It  has  been  found  that  the 
location  of  the  poles  must  be  known  to  7  significant  digits  [16].  This  accuracy 
is  easily  obtained,  however,  using  the  interval-halving  method. 


3.4.1  Extraction  of  Surface  Wave  Poles 


As  pointed  out  in  Chapter  2,  surface  wave  pole(s)  can  occur  in  the 
integrands.  It  is  possible  to  evaluate  the  integrals  as  given  in  Chapter  2.  In 
this  case,  8  should  be  on  the  order  of  ,01&o-  However,  evaluating  the  integral 
in  the  neighborhood  of  these  poles  is  very  difficult.  Therefore,  it  is  desirable 
to  develop  a  technique  by  which  these  poles  may  be  extracted  leaving  a  well- 
behaved  integrand.  This  may  be  done  by  taking  the  -1  term  of  a  Laurent  series 
for  the  integrand  about  the  surface  wave  pole.  Also,  as  pointed  out  in  Chapter 
2,  all  of  the  surface  wave  poles  occur  in  the  range  k  <  p  <  tRk.  Therefore  we 
may  write 


00  2,r  \ 

7  [  N{9,p ) 

J  J  Dip) 

o  o  vr/ 


k  2ir 

dpdO  =  J  J 

o  o 


N(0,p) 

D(P) 


dp  dO 
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sftHk  2  r 


N{0,  p)  _  Rtmo(ptmo) 

D{p)  P  -  PTMO 


dp  dd 


*  Jl% 


(3.28) 


*0  -6^ 

0  0 

PTM 0-6  2jt  / 


k0  0 
2j r 


/>TMO+^  0 
oo  2ir 


JS“-(p,0)z„(p,t)J:?ra(e,e)-  R™’"  )  *<<//? 

P  ~  PTMO/ 


Jsu"(p,»)z„(p,«)J:rM  -  ■  R™;‘  ■)  ^ 

p-pTMO/ 


+  (  I  JT\p,  9)Zz*(p,  0)JT^  °)  d«  W 

V^k0  0 

\  V  *0  “  PTMO  J  J 


(3.29) 


k0  0 
V^k0  2  ir 


Jlrtm(p,O)Zzx(p,e)Jl^0)  -  -  dads 

P  —  PTMO  J 


k o  zn 

A-«(ui  =  j  j  3T'(p,f)z„MjT^9)dadl) 

0  0 

PTMO-&  1-K  ,  \ 

+  /  /  •£"•(*  #)z=(/>.  »)  -  r  *> 

4  '  V  p-ptm o) 

V*Rk 0  2jt  ,  \ 

+  11  (jrm{p>o)z«MJTM  ~  do 

P-PTMO/ 

PTMO+fi  u 
oo  2*r 

+  /  / j%*{p,e)Z„{p,t)j%~(p,»)daW 


\/*Rk0  o 

+  RTMl 

\  \  *0  -  PTMO  )  j 


(3.30) 
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ko  2ir 


K,4U)  =  j  j  j<xr-(p, B)z,,(p,0)j'’r,‘(p,e)dadi) 


0  0 
PTMO-S  2 ir 


+ 


+ 


/  /  [J^(p,o)zX2(p,e)r2fpa(p,d) 


k0  0 
s/tRk  0  2ir 


J  J  UTm(p,o)zxz(P,e)Jezra(p,o) 


ptmo+6  0 
oo  2r 


RtMxz  j 
p  ~  Ptmq) 

RtMxz  \ 
P  —  Ptmq) 


da  df3 

da  d(3 


+  /  l }T'(p,«)z,.MjTM^dl3 

v/JSfal  o 

+  fin  f  ^  ~  -  i*)  Rtu,. 

\  \  k0  —  PTMQ  J  J 


(3.31) 


kr,  2  n 


=  f  J  jT-(p,e)z„(p,e)j~’‘(P,e)dadi3 


0  0 
PTMO-S  2ir 


p  ~  Ptmq 

Rtmxx 
P  —  PTMQ, 


'M0-»  i*  /  p  \ 

/  /  (jy-(p, °)Z**(P,  e)Jl?\p,  0)  -  ™"  ' 

fco  0  ' 

s/?Rko  2ir  , 

I  /( 

PTMO+®  0 

oo  2  ir 

+  /  j  3T'MZ^)j7\P^)dadl3 


+ 


+ 


JT’(p^)z^p^u:7a(p^)- 


da  d/3 

da  d3 


\/<RkQ  0 


+ 


(in  (^o-PTMo\t\  Rj 

V  \  k0-  ptmq  /  / 


Mix 


(3.32) 


where 


2t 

Rtm„  =  j  9)Z„tm(p,  B)j\Y\p, »)  dO 

0 

2ir 

Rtmzx  =  J  JT'iP,  0)ZzxTM{p,  0)Jexx/a(p,  9)  de 
0 
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2ir 

Rtu„  =  J 
0 

2ir 

«™»  =  /  3%*(p,t)Z*rM{p,i)%r(p,f)M. 

0 

(3.33) 


3.5  Extraction  of  asymptotic  forms 

The  second  difficulty  encountered  in  the  evaluation  of  the  integrals  is 
the  evaluation  of  the  asymptotic  portion  of  the  integral.  The  integrand  decays 
slowly,  as  will  be  shown,  and  therefore  the  integration  must  be  carried  out 
to  a  very  large  value.  This,  too,  is  numerically  expensive.  Therefore,  it  is 
desirable  to  extract  the  asymptotic  portion  of  the  integrand  and  integrate  it 
analytically.  However,  this  has  only  been  satisfactorily  implemented  for  the 
one-dimensional  case.  An  alternative  method  has  been  implemented  for  the 
two-dimensional  case. 

3.5.1  One-dimensional  case 

In  the  calculation  of  the  propagation  constant  and  the  current  distri¬ 
bution  for  the  modes  of  a  uniform  line,  it  is  possible  to  obtain  an  asymptotic 
form  of  the  integrand,  which  may  be  integrated  in  closed  form.  We  now  exam¬ 
ine  the  asymptotic  form  of  the  Green’s  function  for  large  a  and  j3.  For  large 
values  of  a,  Equations  2.37-2.50  become 

7i  ~ 


72  ~ 


(3.34) 

(3.35) 
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YtM\ 

YtM2 


jue0eR 


jue0eR 


jue0(  1  +  cr)’ 

jui* 

2a  ’ 

and 

a 


(3.36) 

(3.37) 

(3.38) 

(3.39) 

(3.40) 

(3.41) 

(3.42) 

(3.43) 

(3.44) 

(3.45) 

(3.46) 

(3.47) 


Therefore,  for  large  values  of  a,  the  spectral  domain  Green’s  function 
components  behave  as  follows. 


_ _ _ _  J^f1  yasymp 

ajue0(  1  +  eR)  2a 

_ _ 2 _ i  —  y^sym-p 

a>o(l+e«)  2aj  “  ’ 

- rf - r  +  ^)  =  Z“fmp,  and 

a  ;we0(l+c/?)  2a 

_ £ _ $  _  rjasymp 

_;’a;eo(l+ffl)  2q3  xx 


(3.48) 


(3.49) 

(3.50) 


(3.51) 
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So  the  Green’s  function  behaves  asymptotically  as 


zzz 

~  0(4), 

a 

(3.52) 

Zzx 

~  0(1), 

(3.53) 

Zxz 

~  0(1),  and 

(3.54) 

Zxx 

~  0(a). 

(3.55) 

Therefore,  in  order  to  ensure  the  convergance  of  the  integrals,  the 
basis  functions  used  to  represent  Jz  must  behave  asymptotically  as  and 
the  basis  function  used  to  represent  Jx  must  behave  asymptotically  as  if 
Galerkin’s  method  is  used. 

The  Fourier  transforms  of  the  Maxwellian  basis  functions  are  given 

by 


Jzm(o)  — 


irw 


Mlrf  +  ~  1)*\) 


,  wa 


+  M\—  ~  (m_  i)*!) 


,  and 


nw 

>(a)  =  17 


T  ..wa  ..  r  ..wa  . 

Ml  —  +  -  Ml—  “  m7rl) 


(3.56) 

(3.57) 


where  Jq(x )  denotes  the  zeroth-order  Bessel  function.  This  function  is  evalu¬ 
ated  numerically  using  a  rational  function  approximation.  From  the  asymptotic 
form  of  Jo  [24],  it  can  be  seen  that  the  asymptotic  forms  for  the  Maxwellian 
functions  are 


J~(a)  ~  -A=sm(^~)  (-1)”-1  =  JZF’fa),  and  (3.58) 
y/TTWa  \  L  / 

JUa)  -  -==—  Sin(^)(-1)”>+'  =iITn"(o).  (3.59) 
y/nwa  wa  \  2  / 
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Therefore,  the  Fourier  transforms  of  the  Maxwellian  functions  behave  asymp¬ 
totically  as 


J im(^) 


J rn  (  ®  ) 


(3.60) 

(3.61) 


The  integrands  then  decay  as  This  is  a  fairly  slow  convergance.  Without 
any  acceleration  technique,  the  integrals  must  be  carried  out  to  a  large  value 
of  a  as  much  as  lOOOfco-  However,  since  the  asymptotic  forms  of  the  basis 
functions  and  Green’s  function  are  relatively  simple,  it  is  possible  to  evaluate 
the  asymptotic  portion  of  the  integral  analytically. 

OO 

a 

-  f)  (tz)  Q  +  Z  ™{aw)  -  wCl{aw])  (3'62) 


j  j^mr(a)Z^’ (a)  j"^(a)  da  =  ((-1)”+") 


0 


^a;eo(l  -f  e/i) 


(3.63) 


/  Jzr,w2zr’(a)Jzr’(a)<ia  =  (<-i  r+") 

a 

(^o(ittR))  O  (;+ -  “Ci(<,K,>)  <3-64> 


uu 

J  jzr,(a)2zr’(a)j£r,(a)da  =  ((-l)”+"«) 

fegTHri)  (^)  (a +  Zsia(aw)  -  wCi(aw) ) 


(3.65) 
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where 


OO 

Ci .(*)  =  -/ 


cos (t) 


(3.66) 


The  cosine  integral  may  be  efficiently  and  accurately  approximated  with  a 
rational  function  [23].  So  the  integrals  are  evaluated  as 


oo  a 

J  Jz(a)Zzz(a)jz(a)da  =  J  Jz(a)Zzz(a)Jz(a)  da 
o  o 

OO 

+  J[Jz(a)Zzz(a)  Jz(a) 


-  jazaymp{a)Zazaymp{a)  J“4ymp(a)]  da 

+  (-i  )«+"-*  f _ 2L - A)  (A) 

'  '  Aeo(l  +  eR)  2  /  \irwJ 

+  —  sin(au;)  —  u>Ci(au;)^  ,  ( 


+  (-1) 


(3.67) 


OO  O 

J  Jz{a)Zzx{a)Jx{a)  da  -  J  Jz(a)Zzx{a)Jx(a)  da 
o  o 

OO 

+  J[Jz(a)Zzx(a)Jx(a) 


-  J°*yrnp(a)Z™ymp{a)J°3ymp(a)}  da 

+  (-i  r+*( — t — VA) 

1  }  \u>eo(l+eR)J  U*/ 

-f  -sin(atu)  —  wCiiaw)^  , 

oo  a 

J  jx(a)Zxz(a)  Jz{a)  da  =  J  jx(a)Zxz(a)jz(a)  da 
o  o 

OO 

+  J[Jx(a)Zxz(a)Jz(a) 


(3.68) 


-  j°ayrnp(a)Z*lymp(a)Jz3ymp(a)]  da 

+  (-i  )*+*( — L — )(%*) 

+  1  \ue0(\  +tR))\w*J 


;)©) 
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OO 

J"  J x ( ^ ) Zj xx doc  — 
0 

+ 


+ 


H — sin(aiy)  —  wC'i(aw ,  and 

a 

J  Jx (^Ot')Zxx{^Ot}J x[ot}  dot 

0 

oo 

/  [Ma)Zxx(a)Jx(a) 

a 

JZ3vmp(a)Z£ymp(a))J?ymp(a)]da 

,  Mm+n+2  (  1  \  /87rmn' 


\ue0{l  +  eR)  J  \  w3 

( —  +  —  sin(atu)  —  tuC^au;)^  . 

\a  a  / 


(3.69) 


(3.70) 


3.5.2  Two-dimensional  case 

In  Section  3.1.2,  it  was  seen  that  the  longitudinal  dependance  of  the 
current  is  represented  by  two  types  of  basis  functions,  piecewise  linear  and 
sinusoidal.  The  piecewise  linear  function  is  repeated  here  for  convenience. 

•ww = 

The  Fourier  transform  of  the  piecewise  linear  function  is 

JpwLn(P)  =  ~  cos(0h))exp(j0zm). 

As  can  be  seen,  the  transforms  behave  asymptotically  for  large  /?  as 

jpWLn(P)  ~  0{—).  (3.71) 

As  noted,  the  sinusoidal  basis  functions  are  truncated.  The  length  of  the  sinu¬ 
soid,  as  pointed  out  by  [15]  should  be  an  integer  number  of  half- wavelengths 
in  order  to  obtain  the  fastest  convergence  of  the  integrals.  This  can  be  seen  as 
follows. 


Let 
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,  .  _  J  sin (0zz)  if  0  <  2  <  zl,  and 
sm\z)  |  q  otherwise. 


Then 


For  zt  — 


lin(0) 


For  large  values  of  0, 


1  1 

2(3  +  02 

1  1 

20-02 

-  5JT^cos((/3  +  A)z,) 

+  \jh !>')*>) 

-  +  A)z,) 


•W/?) 


1 

1 

2  0 

+  0Z 

1 

1 

2/? 

-A 

1 

1 

20 

+  0Z 

1 

1 

20 

-A 

1 

1 

2  0 

+  0Z 

1 

1 

2/? 

-0z 

cos(/?z<) 

COs(/?2j) 

sin(0zt) 
sin  (0zt). 


(3.72) 


(3.73) 


(3.74) 


r(l  +  cos  (0zt)  +  sin  (0zt)). 


(3.75) 
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Therefore, 

l,n  ~  O(i).  (3.76) 

On  the  other  hand,  if  the  sine  function  were  truncated  at  an  odd  integer  number 
of  quarter-wavelengths,  FJtn(/?)  would  decay  as  j  for  large  /?.  It  is  desirable  to 
truncate  the  sine  function  at  an  integer  number  of  half-wavelengths  not  only 
because  this  speeds  convergence,  but  also  because  the  subsectional  and  entire 
domain  basis  functions  will  then  decay  at  the  same  rate.  As  it  stands,  J3,„(5) 
can  have  a  pole  at  0  =  0Z.  However,  by  further  constraining  the  truncation  to 
an  even  integer  number  of  half  guide  wavelengths,  we  can  eliminate  the  pole. 
Let 


zl  — 


(3.77) 


Then 


Therefore, 


cos (/3zl,)  —  1,  and 
sin  {/3zl)  =  0. 


J„n(3)  =  0 


(3.78) 


(3.79) 


As  noted  in  Chapter  2,  the  traveling-wave  function  does  not  have  to 
be  truncated  at  all  in  order  for  its  Fourier  transform  to  exist.  That  is 


J source  (^)  — 


e  if  z  <  0,  and 

0  otherwise. 


(3.80) 
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The  Fourier  transform  is  then  derived  as  follows. 

o 

lourM  =  /  e~i0'*e’0*  dz 

— oo 

=  J  ei(P-P*)z  dz 


Therefore, 


—  OO 
0 


=  J  cos(/3  —  [3z)zdz  +  j  J  sin(/?  —  /3z)z  dz 

-oo  -oo 

.  +oo  0 

=  2  /  +  j  /  sin(/?  —  0z)z  dz 


— OO 

+oo 


=  ^  J  ^2dz  +  jlim  J  e“2sin (ft  —  fiz)z  dz 


1 


Ji0urce(^)=2[2^-^)] 


+j  lim 

a — >0 


1  o 


L«2  +  (/?  ~  A)2 


(asin([/?  —  ftz)z)  -  (ft  -  ftz)  cos([ft  -  ftz]z)) 


lource(ft)  =  *6(13  -  ftz)  +  ^  ■  (3.81) 

It  is  immediately  seen  that  this  function  presents  two  numerical  problems. 
First,  it  has  a  pole  at  ft  =  ftz.  Also,  it  behaves  asymptotically  as  Therefore, 
it  is  numerically  more  expensive  to  evaluate  than  the  source  function  used 
here.  However,  it  should  be  noted  that  it  is  more  rigorously  correct  and  may 
be  necessary  when  surface  wave  coupling  is  strong  or  complex  circuit  geometries 
generate  strong  coupling. 

For  large  values  of  ft,  the  spectral  domain  Green’s  function  compo- 
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nents  behave  as  follows. 


So  the  Green’s  function 


/3 _ a2jtjJH 

jue0(\  +  6r)  2a3 

a,  P  jup, 

PKjut0{  1+cr)  2/3  h 

a,  P  ,  , 

~  /3(jW1  +  £„)  +  2/3 

_ _ jufi 

/3jue0(l  +  £r)  2/3 

behaves  asymptotically  as 


(3.82) 

(3.83) 

(3.84) 

(3.85) 


Zzz 

~  0(13), 

(3.86) 

Zzx 

~  0(1), 

(3.87) 

Zxz 

~  0(1),  and 

(3.88) 

Zxx 

(3.89) 

From  the  asymptotic  forms  of  the  Green’s  function  for  large  f3  and  the  asymp¬ 
totic  forms  of  the  basis  functions,  we  see  that  the  two-dimensional  Zzz  integrals 
converge  as  the  Zzx  and  Zxz  integrals  converge  as  and  the  Zxx  integrals 
converge  as  -p-  along  the  j3  axis. 

In  the  two-dimensional  case,  the  extraction  of  the  asymptotic  form  is 
not  straightforward.  Because  the  transverse  dependence  of  the  two-dimensional 
basis  functions  is  the  same  as  that  of  the  one-dimensional  basis  functions  used  in 
the  analysis  of  uniform  line,  the  convergence  along  the  a  axis  is  the  same  as  that 
for  the  one  dimensional  analysis.  Along  the  /3  axis,  the  convergence  is  better 
because  the  basis  functions  exhibit  more  rapid  convergence.  The  integration 
is  performed  in  a  cylindrical  coordinate  system,  and  while  an  asymptotic  form 
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is  available  for  the  Green’s  function  for  large  p,  no  such  forms  are  available 
for  the  basis  functions.  This  is  because  they  are  separable  in  a  and  /?.  The 
asymptotic  form  of  the  Green’s  function  for  large  p 1  is  given  by 


Zzz(p,0) 


Zzx{p,  0) 


Zxz(p,0) 


Zxx(Pi 0) 


pcos2(0) 

ju>p  sin2(0) 

(3.90) 

ju>€o(  1  +  eR) 

2  P 

132  t 

jujpa2 

(3.91) 

jue0(l  +  eR)p 

2  p3  ’ 

cos(0 ) sin(fl) 

juip  cos (9)  sin (9) 

(3.92) 

ju>e0(l  +  eR)p 

2pz 

cos (6)  sin(0) 
jue0(l  +  tR)p 

jujpcos(d)  sin(0) 

2„3  ’  and 

(3.93) 

psin2(0) 

ju>p  cos2(0) 

(3.94) 

jueQ(l  +  eR) 

2  P  ' 

For  a  conductor  backed  homogeneous  medium  of  relative  dielectric 
constant  ee,  the  space  domain  Green’s  function  is  given  in  closed  form  by 


Zzz(x  -  x0,z  -  z0 )  = 

-Mo  +  £r)[ 


~J  f.  1.2  ,  &  ueXP (~JkRz) 


ivueo  tR 


dz2‘ 


R, 


exp  (—jkRi) 

% 


],  (3.95) 


Zzx(x  Xq,  Z  Zq )  — 

-j 


2  ,  d2  Nr  exp  i-jkR,) 


j  (  .  L2  ,  ^  J' 

\irujtQtR  e  0  dzdx  R, 


exp  (-jkRj) 

Ri 


],  (3.96) 


ZXz{x  z  zq)  — 

-(a°  +  a^)[ 


-J  /_  ,.2  ,  &  \  r  exp(  j  ^  R$ ) 


47ru;eoeK 


R, 


exp  j-jkRt) 

Ri 


],  and  (3.97) 


^he  asymptotic  form  for  Zzz  in  [36]  appears  to  be  incorrect. 
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*EQ  y  Z  Zo)  ■“ 

I 


J  /.  ,.2  ,  d2  xrexP(-i^»)  exp (-jkRi) 


4irue0eR 


R, 


Ri 


where 


Ra  =  \J(z  -  20)2  +  (x  -  x0)2,  and 
Ri  =  \J(z  -  z0)2  +  (x  -  x0)2  +  (2d)2. 


(3.98) 

(3.99) 
(3.100) 


Using  the  identity  given  on  page  179  of  [25]  and  the  differentiation  theorem  for 
Fourier  transforms,  the  spectral  domain  form  of  this  Green’s  function  is  given 
by 


Zzz(a,0)  = 

2]ue0eey2 

(3.101) 

Zzx{a,0 )  = 

(2;W,7j>(1  exp(  272 d)) ’ 

(3.102) 

II 

2,W,7!)(1  eXP(  27jrf))'and 

(3.103) 

ZXT(a,0)  = 

2jwe  ote72 

(3.104) 

This  behaves  asymptotically  for  large  p  as 

Z„(a,0)  =  $- 


S' 


Zzx(a,/3)  = 


2  p  2ju>tottP' 

kl  0a 


—  o 


2  p  2juc0eep' 


k2 


Zx,(a,(3)  =  r*- 


aS 


2  p  2jut0eep 

k2  a2 


,  and 


Z«(a,0)  =  r^-rr 


2  p  2jue0eep 

Therefore,  if  we  let  ce  =  the  asymptotic  forms  become 


(3.105) 

(3.106) 

(3.107) 

(3.108) 


Zzz{oc,0)  = 


jup 


jue0(l+eR)p  2  p' 


(3.109) 
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! 3a 

jup 

ju}e0(l  +  eR)p 

2  P  ’ 

a/3 

-  “d 
2  P 

jue0(l  +  tR)p 

a2 

jup 

jut  0(1  +  tR)p 

2  P  ‘ 

(3.110) 

(3.111) 

(3.112) 


It  can  be  seen  that  the  TM  term  of  the  asymptotic  form  for  the  homogeneous 
case  is  identical  to  that  of  the  grounded  dielectric  slab.  Therefore,  this  aymp- 
totic  form  may  be  subtracted  from  the  spectral  domain  Green’s  function  to 
yield  a  more  rapidly  converging  kernel.  Although  the  additive  term  may  not 
be  evaluated  analytically  as  in  the  one-dimensional  case,  it  may  be  evaluated  in 
the  space  domain.  This  requires  the  computation  of  a  four-dimensional  finite 
integral.  When  simple  basis  functions  are  used,  this  four-dimensional  integral 
may  be  further  reduced  to  two-dimensions  because  two  of  the  integrations  may 
be  performed  analytically.  However,  the  use  of  the  Maxwellian  functions  allows 
only  one  of  the  integrations  to  be  done  analytically.  Because  the  integration 
is  in  the  space  domain,  a  singularity  exists  at  the  source  point.  In  order  to 
perform  the  integration  efficiently,  this  singularity  must  be  extracted  and  inte¬ 
grated  analytically. 


s. 

OO  2 


/  /  J,r,-{p,»)z„(p,e)jr‘(p,0)pdpM  = 


0  0 


2L 

OO  2 


/  J  9)  [Z„(/>,  0)  -  Zk„(p,  0)J  j\^(p,  0)pipd8  + 


0  0 


an 

X  Z  XQ  ZQ 


Jt2at(x,z)Zzz(x  —  xq ,  2  —  zo)Je2xpa(x,  z)dxodzodx  dz(3.\]3) 


Because  the  space  domain  integration  is  still  somewhat  laborious,  it  nas  been 
found  that  this  method  did  not  reduce  the  computation  time  significantly. 
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3.6  Numerical  Evaluation  of  the  Integrals 

Because  the  the  range  of  integration  is  the  same  for  all  the  integrals, 
it  is  advantageous  to  precompute  and  store  the  components  of  the  Green’s 
function  and  the  expansion  and  testing  functions.  Thus  the  Green’s  function 
is  computed  only  once  for  each  integration  step  for  the  entire  matrix.  This  has 
been  noted  by  [2].  Some  savings  is  realized  in  the  computation  of  the  basis 
functions  due  the  eliminated  redundancy.  This  is  particularly  suited  to  vector 
processors  because  the  computation  of  the  integrals  may  be  vectorized. 


Chapter  4 


THE  OPEN-CIRCUIT  MICROSTRIP 
DISCONTINUITY 

4.1  Introduction 

Microstrip  open-end  discontinuities  are  an  integral  part  of  most  mi¬ 
crostrip  circuits  and  antennas.  The  geometry  of  a  microstrip  open-end  discon¬ 
tinuity  is  shown  in  Figure  4.1.  Open-end  discontinuities  in  open  microstrip  are 
imperfect  open-circuits  because  of  energy  storage  in  non-radiating  fields  and 
energy  leakage  into  radiating  fields.  Energy  leakage  into  both  space  waves  and 
surface  waves  occurs.  The  discontinuity  exhibits  a  minimum  susceptance  immi- 
tance  function.  Therefore,  an  appropriate  circuit  model  is  a  parallel  combina¬ 
tion  of  a  reactance  and  a  conductance  as  shown  in  Figure  4.2.  In  general,  both 
the  reactance  and  the  conductance  are  frequency  dependent.  In  this  chapter, 
microstrip  open-end  discontinuities  in  open  environments  are  analyzed  using  a 
combination  of  the  the  deterministic  and  eigenvalue  formulations  of  the  spec¬ 
tral  domain  approach  presented  in  Chapters  2  and  3.  The  study  includes  the 
effects  of  both  longitudinal  and  transverse  currents  and  incorporates  the  proper 
edge  conditions  for  each  component  of  the  current  at  the  edges  of  the  strip  and 
the  open  end  of  the  strip. 
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groundplane 


Figure  4.1:  Geometry  of  Microstrip  Open-End  Discontinuity 


z=0 


Figure  4.2:  Equivalent  Circuit  of  Microstrip  Open-End  Discontinuity 
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4.2  Method 

The  computational  method  consists  of  two  parts.  In  the  first  part, 
the  spectral  domain  eigenvalue  formulation  is  used  to  compute  the  propagation 
constant,  characteristic  impedance,  and  longitudinal  and  transverse  current 
distributions  for  the  fundamental  mode  on  an  infinitely  long  open  microstrip 
corresponding  to  the  uniform  region.  The  source  current  to  be  used  in  the 
second  part  of  the  method  is  then  given  by 

L 

=  £«/•/*/(*)  exp(-j&.z),  and  (4.1) 

i=i 
L 

=  ^2btJxi(x)exp{-jP2z).  (4.2) 

/=i 

The  functions  Jzm  and  Jxn  are  the  Maxwellian  basis  functions  described  in 
Chapter  3  They  include  the  proper  edge  condition  for  the  longitudinal  and 
transverse  current  at  the  edges  of  the  microstrip.  The  constants  ai...ai  , 
bx . . .  b £,  and  /?,  the  propagation  constant,  are  determined  in  the  spectral  do¬ 
main  analysis.  In  practice,  L  —  2  yielded  good  results  and  was  therefore  used 
throughout  the  analysis.  The  propagation  constant  for  the  fundamental  mode, 
(3Z,  is  always  greater  than  the  value  of  p  for  any  surface- wave  pole.  Therefore, 
the  fundamental  mode  is  always  slow  compared  with  any  surface  wave  mode 
and  thus  does  not  radiate  any  power  into  surface-waves.  In  the  second  part  of 
the  method,  the  microstrip  is  divided  into  two  overlapping  regions:  a  uniform 
region  far  away  from  the  open  end  and  a  perturbed  region  near  the  open  end, 
as  in  [15].  In  the  uniform  region,  the  current  is  assumed  to  consist  of  only  the 
fundamental  mode.  This  will  be  true  if  the  higher  order  modes  are  either  leaky 
or  evanescent.  The  current  in  the  uniform  region  is  represented  by  the  source 


jrrce(*,z) 

JTrce(*,z) 
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function  and  a  reflected  wave  with  an  unknown  amplitude,  T. 

L 

J“n‘/°rm(x,  z)  =  ^2aiJzi(x)(exp(—j 0zz)  —  T  exp(+j/?2z)),  and  (4.3) 

1=1 

=  'tblJxl(x)(exV(-jPzz)  +  rexp(+jpzz)).  (4.4) 

i=i 

Note  that  T  is  the  voltage  reflection  coefficient.  As  discussed  in  Chapter  3, 
the  exponential  functions  are  split  into  sine  and  cosine  functions  which  are 
truncated  to  give  the  most  rapidly  converging  Fourier  transforms. 

The  longitudinal  dependence  of  the  current  in  the  perturbed  region 
near  the  open  end  is  augmented  with  piecewise  linear  (rooftop)  subsectional 
basis  functions,  which  are  defined  in  Chapter  3.  The  piecewise  linear  func¬ 
tions  allow  a  reasonable  approximation  to  the  edge  condition  for  the  both  the 
transverse  and  the  longitudinal  currents  at  the  open  end.  The  transverse  de- 
pendance  of  the  current  in  the  perturbed  region  is  represented  with  Maxwellian 
basis  functions  with  variable  coefficients.  This  takes  into  account  mode  con¬ 
version  in  the  perturbed  region.  Therefore,  the  current  in  the  perturbed  region 
is  given  by 

M  N 

=  E  J*rn(x)  E  CmnFn(z)  +  *),  and 

m=l  n=l 

M  N 

Jr'*’ W(X,  z)  =  £  J,m(x)  £  dmnFn(z)  +  z). 

m= 1  n=  1 

where  F\(z) . . .  Fn(z)  are  piecewise  linear  subsectional  basis  functions.  The 
positioning  of  the  basis  functions  is  shown  in  Figure  4.3. 

If  just  the  longitudinal  current,  Jz,  is  taken  into  account  and  mode 
conversion  is  neglected  ( M  =  1),  as  it  was  in  [15],  there  are  N  +  1  unknowns, 
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Figure  4.3:  Positioning  of  Longitudinal  Basis  Functions  for  Analysis  of  Mi¬ 
crostrip  Open-end  Discontinuity 
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the  complex  coefficients  of  the  subdomain  expansion  functions  and  the  complex 
reflection  coefficient.  Since  it  is  not  possible  to  test  the  equation  with  the  source 
function,  Galerkin’s  method  cannot  be  used  here.  However,  it  is  still  useful 
to  use  the  M  subsectional  basis  functions  as  testing  functions  in  a  “nearly 
Galerkin”  method.  This  generates  an  underdetermined  linear  system  with  one 
less  equation  than  unknowns.  One  method  of  solving  the  problem  is  to  simply 
use  one  more  testing  function  to  generate  one  more  equation.  The  linear  system 
is  then  square  and  may  be  solved  using  Gaussian  elimination.  This  method  was 
proposed  by  [15].  Thus,  the  following  linear  system  is  derived. 


' 

■  J^zzpp 

rszzpp 

Kl,2 

TfZZpp  -| 

...  KlN 

■  Kzzp-  ■ 

■ 

(Zi 

-  _szzp+  ' 
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...  i\2,s 
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ry-zzpp 
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l/zzpp 
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kn 

as 
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‘-’N 

js*zpp 

L  AAT+1,1 

J/*zpp 

nN+l,2 

j/zzpp 

•  •  •  nN+l,N 

TSZZP- 

L  ajv+ i  j 

• 

&N+ 1 

czzp+ 

.  J 

The  portion  of  the  matrix  generated  by  testing  and  expanding  with  subdomain 
expansion  functions  has  Toeplitz  symmetry.  Therefore,  only  one  row  of  ele¬ 
ments  needs  to  be  calculated.  To  generate  the  entire  linear  system,  4Ar  +  2 
integrals  must  be  computed. 

When  two  dimensional  current  is  considered,  the  linear  system  be¬ 
comes 


■  [Kzz]  [Kxx\  ' 

Ml 

f  [S,] 

.  [Kxz]  [Kxx\ 

[b]  j 

[  [5.1 

where  Kzz  denotes  the  coefficient  matrix  given  above  and  Kzx  etc.  are  similar. 
By  reciprocity,  Kzx  =  Kxz.  We  note  that  a^+i  =  —  T  and  bs+i  =  T.  Therefore, 
the  system  of  equations  given  above  is  over-determined  by  one  equation.  The 
overdetermined  system  was  solved  using  a  least-squares  approach.  The  numer- 
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ical  cost  of  including  the  effects  of  two-dimensional  current  are  substantial;  we 
must  evaluate  UN  +  5  two-dimensional  integrals. 

One  drawback  to  this  method  is  that  there  are  several  unknown  quan¬ 
tities  which  must  be  determined  empirically:  the  length  of  the  uniform  line  basis 
functions,  z/,  the  length  of  the  region  of  subdomain  expansion,  and  the  length 
of  the  subdomain  expansion  functions.  It  has  been  empirically  determined 
that  using  z\  —  6Ae//  yields  satisfactory  results.  Also,  it  has  been  found  that 
using  subdomain  expansion  functions  of  width,  2 h  =  -fff,  yields  satisfactory 
results.  However,  the  length  of  the  subdomain  expansion  region,  and  hence 
the  number  of  subdomain  expansion  functions,  depends  on  the  geometry  of  the 
discontinuity. 

Some  of  this  uncertainty  may  be  removed  by  generating  an  overde¬ 
termined  system  by  using  more  subsectional  testing  functions  to  generate  more 
equations.  The  over  determined  system  is  solved  via  a  least  squares  method.  T 
is  the  complex  reflection  coefficient  for  the  fundamental  mode.  From  a\ . .  .a\j 
and  T,  the  current  over  the  entire  strip  can  be  determined.  From  T,  the  reflec¬ 
tion  coefficient,  and  Zo,  the  characteristic  impedance  of  the  uniform  line,  an 
equivalent  admittance  can  be  calculated 

=  Vo  (j^y)  (4.5) 

where  Yq  =  ^r-. 

4.3  Numerical  Results 

In  Figures  4.4  and  4.5,  the  longitudinal  (z-directed)  current  and  the 
transverse  (x-directed)  current  on  an  open  circuited  microstrip  line  are  plotted 
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Figure  4.4:  Magnitude  of  z-directed  Current  on  Microstrip  Open-end  Discon¬ 
tinuity 

to  show  the  edge  conditions  for  a  typical  case.  The  characteristics  are  as  follows: 

•  Relative  dielectric  constant,  cr  :  12.8 

•  Substrate  height,  h  =  .300  mm 

•  Strip  width,  w  =  .600  mm 

•  Frequency  =  4.00  GHz 

It  can  be  seen  that  the  piecewise  linear  functions  allow  the  edge  conditions  at 
the  open  end  for  both  the  transverse  current  and  the  longitudinal  current  to 
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Figure  4.5:  Magnitude  of  x-directed  Current  on  Microstrip  Open-end  Discon¬ 
tinuity 
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open  microstrip  tine: 


10.00  20.00  30.00  40.00 


Figure  4.6:  Effective  Dielectric  Constant  of  Uniform  Microstrip  Line 
be  accurately  modeled. 

In  order  to  allow  comparison  with  previously  published  experiment  al 
data,  [49]  an  analysis  was  made  of  an  open-end  discontinuity  in  microstrip  line 
with  the  following  parameters. 

•  Relative  dielectric  constant,  :  9.9 

•  Substrate  height,  h  =  .635  mm 

•  Strip  width,  w  =  .600  mm 

The  dispersion  curve  for  ce//  is  shown  in  Figure  4.6.  In  Figure  4.7,  the  char¬ 
acteristic  impedance  is  plotted  as  a  function  of  frequency.  Figures  4.8  and  1.9. 
show  the  magnitude  and  phase  of  the  reflection  coefflcent  for  an  open  circuit 
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Figure  4.7:  Characteristic  Impedance  of  Uniform  Microstrip  Line 

discontinuity.  In  Figures  4.10  and  4.11  G  and  B  are  plotted  as  functions  of 
frequency.  This  line  is  very  narrow  electrically;  at  20  GHz  the  line  is  still 
only  0.04  free-space  wavelengths  wide.  The  phase  agrees  reasonably  well  with 
the  measured  data  in  [49]  (the  magnitude  was  not  presented  in  [49])  and  with 
calculated  data  in  [18]  and  [44]. 

In  order  to  quantify  the  efFect  of  the  transverse  current  on  the  reflec¬ 
tion  coefficient  of  a  microstrip  open-end  discontinuity,  an  analysis  was  made  of 
an  open-end  discontinuity  in  microstrip  line  with  the  following  parameters. 

•  Relative  dielectric  constant,  ir  :  9.6 

•  Substrate  height,  h  =  .635  mm 

•  Strip  widths,  w  —  .600  mm  and  w  =  1.20  mm. 
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Figure  4.10:  Equivalent  Conductance  of  Microstrip  Open-end  Discontinuity 
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Figure  1.11:  Equivalent  < 'apacitaru  e  of  M irrost  rip  Open  end  Discontinuity 
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Effective  Dielectric  Constant 


Figure  4.12:  Effective  Dielectric  Constant  of  Uniform  Microstrip  Line 

The  narrow  line  is  very  similar  to  that  in  the  above  analysis.  It  was  chosen 
because  it  is  a  very  commonly  used  design.  The  dispersion  curves  for  ee// 
are  shown  in  Figure  4.12.  In  Figure  4.13,  the  characteristic  impedances  are 
plotted  as  a  function  of  frequency.  As  can  be  seen,  the  narrow  line  exhibits  a 
characteristic  impedance  which  varies  from  approximately  50  ohms  at  10  GHz 
to  70  ohms  at  40  GHz.  The  characteristic  impedance  of  the  wide  line  varies 
from  approximately  35  ohms  to  45  ohms  over  the  same  range.  In  Figure  4.14. 
the  ratio  of  the  transverse  to  longitudinal  current  is  plotted.  It  can  be  seen  that 
the  transverse  current  on  the  narrow  line  is  not  very  significant  in  comparison 
to  the  longitudinal  current.  However,  the  transverse  current  on  the  wide  line 
can  be  quite  appreciable.  In  Figure  4.15  and  Figure  4.16  the  magnitude  and 
the  phase  of  the  reflection  coefficient  of  the  open-end  discontinuity  are  plotted. 
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Figure  4.15:  Magnitude  of  Reflection  Coefficient  of  Microstrip  Open-End  Dis¬ 
continuity  j  =«  .945 

It  can  be  seen  that  the  inclusion  of  the  transverse  current  in  the  analysis  mainly 
affects  the  phase  of  the  reflection  coefficient.  It  can  also  be  seen  that  the  effect 
of  the  transverse  current  is  a  decrease  in  the  phase  of  the  reflection  coefficient. 
This  seems  resonable  because  the  transverse  current  tends  to  store  magnetic 
energy.  The  longitudinal  current  stores  mainly  electric  energy.  Thus  the  effect 
of  the  transverse  current  can  cancel  the  effect  of  the  longitudinal  current  to 
some  extent. 

In  order  to  further  compare  this  analysis  with  previous  numerical 
work,  and  analysis  was  made  of  an  open-end  discontinuity  in  microstrip  line 
of  width  .lAo  on  a  substrate  with  a  relative  dielectric  constant  of  12. S.  The 
magnitude  of  the  reflection  coefficient  of  the  open-circuited  microstrip  line  is 
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Figure  4.17:  Magnitude  of  Reflection  Coefficient  of  Microstrip  Open-End  Dis¬ 
continuity  ^  ss  1.90 
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Figure  4.19:  Magnitude  Reflection  Coefficient  Coefficient  of  Microstrip  Open- 
end  Discontinuity 

shown  in  Figure  4.19.  In  the  same  figure,  the  magnitude  of  the  reflection 
coefficient  is  shown  calculated  including  only  longitudinal  current  on  the  strip 
as  was  done  in  [15].  This  data  agrees  well  with  that  presented  in  [15].  It  can 
be  seen  that  for  narrow  lines  such  as  this  one,  the  inclusion  of  the  transverse 
current  makes  little  difference  in  the  magnitude  of  the  reflection  coefficient  for 
small  values  of  However,  near  the  cutoff  of  the  first  TE  surface  wave,  it  is 
seen  that  the  magnitude  of  the  reflection  coefficient  displays  a  more  pronounced 
trough  and  peak.  This  seems  reasonable  because  the  inclusion  of  the  energy 
stored  in  the  transverse  current  could  increase  the  frequency  sensitivity  of  the 
reflection  coefficient. 
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Also  in  Figure  4.19,  the  magnitude  of  the  reflection  coefficient  for 
an  open-circuited  microstrip  line  of  width  .15Ao  on  a  substrate  with  a  relative 
dielectric  constant  of  12.8  is  shown.  It  can  be  seen  that  the  wider  line  has  a 
more  pronounced  resonance  at  the  cuton  of  the  first  TE  surface  wave  mode. 

4.4  Conclusions 

Microstrip  open-circuit  terminations  have  been  analyzed  using  a  de¬ 
terministic  spectral  domain  method.  The  radiation  loss  from  the  open-end 
discontinuity  can  be  significant  when  the  substrate  is  electrically  thick.  It  has 
been  shown  that  the  the  transverse  current  in  the  microstrip  has  little  effect 
on  the  reflection  coefficient  as  long  as  the  line  is  narrow  and  the  substrate  is 
electrically  thin.  However,  the  effects  of  the  transverse  current  are  more  pro¬ 
nounced  for  wider  lines  and  thicker  substrates.  The  effect  is  most  noticeable 
near  the  cutoff  frequencies  of  the  TE\  surface- wave  mode. 


Chapter  5 


THE  SHORT-CIRCUIT  SLOTLINE 
DISCONTINUITY 


5.1  Introduction 

The  slotline  short-circuit  discontinuity  is  the  simplest  slotline  discon¬ 
tinuity.  The  geometry  of  the  slot  line  short-circuit  is  shown  in  Figure  5.1.  Like 
the  open-circuit  discontinuity  in  microstrip,  the  short-circuit  slotline  disconti¬ 
nuity  is  non-ideal  because  of  energy  storage  near  the  discontinuity  and  energy 
loss  from  radiation  into  surface  and  space  waves.  The  slotline  short-circuit 
discontinuity  exhibits  a  minimum  reactance  immitance  function.  Therefore,  an 
appropriate  equivalent  circuit  for  a  short-circuit  slotline  discontinuity  is  the  se¬ 
ries  combination  of  a  resistance  and  an  inductance  both  of  which  are  in  general 
frequency  dependent  as  shown  in  Figure  5.2. 

In  this  chapter,  slotline  short-circuit  discontinuities  in  open  environ¬ 
ments  are  analyzed  using  an  extension  of  the  method  presented  in  Chapters 
2,3,  and  4.  This  study  includes  the  effects  of  both  longitudinal  and  transverse 
electric  slot  fields  and  incorporates  the  proper  edge  conditions  for  each  com¬ 
ponent  of  the  field  at  the  edges  of  the  slot  and  the  short-circuited  end  of  the 
slot. 

Although  the  microstrip  open-circuit  and  the  slotline  short-circuit 
are  similar  in  that  both  radiate  power  into  space  waves  and  surface  waves,  the 
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Figure  5.1:  The  Slotline  Short-Circuit  Discontinuity 


z=0 


Figure  5.2:  Equivalent  Circuit  of  Slotline  Short-Circuit  Discontinuity 
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excitation  mechanisms  are  different  and  therefore,  the  radiation  characteristics 
are  different.  In  this  chapter  we  compare  the  non-ideal  aspects  of  open-circuit 
microstrip  and  short-circuit  slotline  discontinuities. 

5.2  Method 

In  principle,  the  method  of  analysis  described  in  Chapters  2  and  3  may 
be  used  to  model  slotline  discontinuities.  However,  because  the  current  exists 
over  semi-infinite  planes,  the  representation  of  the  unknown  current  with  basis 
functions  would  be  difficult  and  inefficient.  Therefore,  an  alternative  approach 
is  used.  Consider  the  spectral  domain  equations  defining  the  Green’s  function. 


Ez(a,0)  =  Zzz(a,0)Jz(a,0)  +  Zzx(a,  0)Jx(a,  0)  (5.1) 

Ex{a,  0)  =  Zxz(a,  0)Jz(ct,  (3)  +  Zsx(a,  0)  J*(a,  0)  (5.2) 

This  expression  can  be  inverted  to  give 

U<*,0)  =  Yxx(a,0)Ex(a,0)  +  Yxz(a,  0)Ez(a,  0)  (5.3) 

•/*(<*,  0)  =  Yzx{a,  0)Ex(a,  0)  +  Yzz(a,  0)Ez{a,  0).  (5.4) 


Because  the  slotline  has  no  ground  plane,  the  Green’s  function  is  slightly  dif¬ 
ferent  from  that  used  for  microstrip  line.  The  Green’s  function  is  derived  in 
Appendix  2  and  is  given  by 


Yr:  = 

-N]Ye  -  N]Yh , 

(5.5) 

Yxz  = 

-NZNX(Y'  -  Yh ), 

(5.6) 

Yzx  = 

-NzNx(Ye  -  Yh)  and, 

(5.7) 

Yzz  = 

V 

1 

t) 

1 

(5-8) 
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where 


Ve  =  Yt  M2 


YtM2  +  Ytm\  h(72«0 ' 


,Ytm\  +  Yjm2  coth(7 2d); 

yh  _  Y  (  ^ TE 2  +  ^ TE1  COtM72^)  \ 

V  Ytei  +  Yte2  coth(72 d) ) 


,  and 


(5.9) 

(5.10) 


Like  the  Green’s  function  used  in  the  analysis  of  microstrip,  this  Green’s  func¬ 
tion  has  poles  corresponding  to  surface  waves  in  the  dielectric.  These  surface 
waves  correspond  to  those  in  a  grounded  dielectric  slab  just  as  those  generated 
by  microstrip  discontinuities.  Therefore,  the  pole  locations  are  the  same  as 
those  for  the  grounded  slab.  However,  the  residues  at  the  poles  are  different 
because  the  surface  wave  excitation  mechanism  is  different.  The  residue  cal¬ 
culation  is  the  same 
replaced  by 


that  given 

in  Chapter  2  except  Zzztm  ■ 

■  •  ZXxT M  are 

YxxTM  = 

(5.11) 

YxzTM  = 

~NxNzYjM 

(5.12) 

Y^xTM  — 

—NxNzYjM  and, 

(5.13) 

YxxTM  = 

~NzYtm, 

(5.14) 

where 


ye 

rTM 


dY{ 

dp 

djl 

dp 


YtM2 

(-P 


/  -  \ 
YtM2  +  YjMl  COth(  72(f) 


V 


3Y'  9Y' 

dp  '  dp 


) 


7i 

P 


,  7i 


2  j  Yfip),  and 

)  yl'(p)- 
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For  TE  modes  the  residue  is  given  by 


YxxTE  =  -N2zY£e, 

(5.15! 

YXzte  =  +NxNzYTkE 

(5.16) 

YxxTE  =  +NxN2Y£e  and, 

(5.17) 

YzzTe  =  ~N2xY£e, 

(5.18) 

where 


Yh 

ITE 


—  YtE2 


Yte2  +  Ytei  coth(72d) 


2Xl.  4. 

dp  '  dp 


(^f)  (*i»  +  rTM2<i72(coth2(7,<()  -  1))  ,  and 

^ ^  +  i'TEMcrth2(^‘‘'>  -  0) 

With  this  formulation  of  the  Green’s  function,  we  can  expand  the  slot  electric 
field  in  terms  of  known  expansion  functions  in  the  same  way  the  strip  current 
was  expanded  for  microstrip.  We  can  make  use  of  the  same  expansion  and 
testing  functions  as  used  in  the  analysis  of  microstrip  except  that  the  functions 
used  to  model  the  transverse  dependance  of  the  current  are  reversed. 


dYf 

dp 

dY2h 

dp 


Exm  ( 2- ) 


Ezm(x) 


cos(2l2Zij»£) 

i/'  -  <£>’ 

sin(^P) 

i/i  -  (V)j 


(5.19) 


(5.20) 


The  functions  Ezm  and  Eim  include  the  proper  edge  condition  for  the  longitudi¬ 
nal  and  transverse  field  at  the  edges  of  the  slot.  The  spectral  domain  eigenvalue 
approach  described  in  Chapter  2  is  used  to  compute  the  propagation  constant 
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and  the  transverse  and  longitudinal  electric  slot  field  components  for  the  fun¬ 
damental  mode  on  an  infinitely  long  open  slotline.  For  the  fundamental  mode, 
the  propagation  constant,  0Z  is  always  greater  than  any  the  value  p  of  any 
surface-wave  pole.  Therefore,  the  fundamental  mode  on  slotline  is  always  slow 
compared  to  any  surface-wave  mode  and  thus  does  not  radiate  into  surface- 
wave  modes.  Once  the  fields  for  the  fundamental  mode  on  uniform  slotline 
are  determined,  the  slotline  is  divided  into  two  regions,  a  uniform  region  far 
from  the  shorted  end  and  a  perturbed  region  near  the  shorted  end  as  in  [15]. 
In  the  uniform  region  the  slot  electric  field  is  assumed  to  consist  of  only  the 
fundamental  mode.  The  slot  electric  field  in  the  uniform  region  is  represented 
by  the  source  function  and  a  reflected  wave  with  an  unknown  amplitude,  I\ 

E™,/orm(x,z)  =  £a,£xl(x)(exp(-j/3zz)  -  V  exp (+jfizz)) 

1=1 

Erform(x,z)  =  Eb,Ezl(x)(exP(-j0zz)  +  Texp(+j0zz)) 

i=i 

The  constants  a\...ai„  bi...bi  and  0,  the  propagation  constant  are  deter¬ 
mined  using  the  spectral  domain  eigenvalue  formulation.  The  longitudinal 
dependance  of  the  slot  electric  field  in  the  perturbed  region  near  the  shorted 
end  is  represented  with  piecewise  linear  (rooftop)  subsectional  basis  functions 
which  are  defined  in  Chapter  3.  The  piecewise  linear  functions  allow  a  reason¬ 
able  approximation  to  the  edge  condition  for  the  both  the  transverse  and  the 
longitudinal  components  of  the  electric  field  at  the  shorted  end.  The  transverse 
dependance  of  the  slot  electric  field  in  the  perturbed  region  is  represented  with 
Maxwellian  basis  functions  with  variable  coefficients.  The  development  of  the 
linear  system  is  then  the  same  as  that  for  the  open-end  microstrip  discontinuity 
which  was  developed  in  Chapter  4. 
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Figure  5.3:  Magnitude  of  Reflection  Coefficient  of  Slotline  Short-circuit  Dis¬ 
continuity 

5.3  Numerical  Results 

In  Figures  5.3  and  5.4,  the  magnitude  and  phase  of  the  reflection  co¬ 
efficient  versus  frequency  is  plotted  for  a  shorted  slotline  on  a  substrate  with  a 
relative  dielectric  constant  of  20.0  and  thickness  of  3.175  mm  for  various  slot 
widths.  In  Figures  5.5  and  5.6,  the  normalized  resistance  and  inductance  of 
the  equivalent  circuit  are  plotted.  It  can  be  seen  that,  as  with  the  microstrip 
open  circuit  discontinuity,  the  radiation  loss  becomes  very  large  near  the  cutoff 
frequency  of  the  second  surface  wave  mode.  In  Figure  5.7,  the  magnitude  of 
the  reflection  coefficient  for  a  short-circuited  slotline  and  an  open-circuited  mi¬ 
crostrip  line  are  plotted  for  comparison.  The  substrate  thickness  and  dielectric 
constant  are  the  same  for  the  microstrip  and  the  slotline.  Also,  the  slot  width  is 


Figure  5.4:  Phase  of  the  Reflection  Coefficient  Slotline  Short-circuit  Disconti 
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R  t  -  20.0 

■ 


Figure  5.5:  Equivalent  Resistance  of  Slotline  Short-circuit  Discontinuity 


x 


Figure  5.6:  Equivalent  Inductance  of  Slotline  Short-circuit  Discontinuity 
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in 


Figure  5.7:  Comparison  of  Slotline  Short-circuit  and  Microstrip  Open-end  Dis¬ 
continuities 
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the  same  as  the  microstrip  width.  The  microstrip  open-circuit  exhibits  greater 
loss  because  the  open-circuited  microstrip  discontinuity  excites  the  T Mq  sur¬ 
face  wave  mode  more  strongly  than  the  slotline  short-circuit  discontinuity.  This 
is  because  the  fields  of  the  fundamental  mode  of  microstrip  match  those  of  the 
T Mq  surface  wave  more  closely  than  those  of  the  fundamental  mode  of  slotline. 


Chapter  6 


THE  MICROSTRIP  GAP  DISCONTINUITY 

6.1  Introduction 

Microstrip  gap  discontinuities  are  commonly  used  to  achieve  capac¬ 
itive  coupling  in  microstrip  circuits  and  antennas.  In  open  structures,  mi¬ 
crostrip  gap  discontinuities  can  radiate  into  space  waves  and  surface  waves 
resulting  in  losses  and  spurious  coupling.  In  this  chapter  we  will  analyze  gap 
discontinuities  and  gap-coupled  resonators.  This  analysis  includes  the  effects 
of  a  two-component,  two-dimensional  current  flow  on  the  conductor  and  mode 
conversion.  Through  the  use  of  the  exact  spectral  domain  Green’s  function 
it  includes  the  effects  of  space- wave  and  surface- wave  radiation.  Like  [4],  this 
method  can  be  applied  to  asymmetric  discontinuities.  However,  it  is  more 
efficient  due  to  the  travelling  wave  source  formulation,  requiring  only  two  ma¬ 
trix  fills  and  inversions  to  obtain  the  two-port  scattering  matrix.  In  the  case 
of  symmetric  structures,  only  one  matrix  fill  and  inversion  is  required  with 
this  method,  whereas  [19]  requires  two.  Unlike  the  discontinuities  analyzed 
in  Chapters  4  and  5,  the  gap  and  gap-coupled  resonator  are  two-port  devices 
and  therefore,  in  addition  to  the  reflection  coefficient  a  transmission  coefficient 
must  also  be  determined. 

In  Figure  6.1  a  microstrip  gap  discontinuity,  gap-coupled  resonator, 
and  coupled  resonator  bandpass  filter  are  shown.  It  was  desired  to  deter- 
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Microstrip  Gap  Discontinuity 
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Microstrip  Gap-Coupled  Resonator 
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Figure  6.1:  Microstrip  Gap  Discontinuity,  Gap-coupled  Resonator,  and  Gap- 
coupled  Bandpass  Filter 
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mine  whether  gap  discontinuities  interact  significantly  on  open  microstrip  when 
spaced  slightly  less  than  one-half  wavelength  apart  as  they  would  be  in  such  a 
filter.  It  is  shown  that  the  gap  discontinuities  do  interact  at  high  frequencies 
reducing  the  Q  and  the  rejection  of  gap-coupled  resonators. 

6.2  Application  of  the  Method  to  Microstrip  Gap  Dis¬ 
continuities 

As  in  the  analysis  of  microstrip  open-end  discontinuities,  the  spectral 
domain  eigenvalue  fomulation  is  used  to  compute  the  propagation  constant 
and  the  longitudinal  and  transverse  currents  for  the  fundamental  mode  on  an 
infinitely  long  open  microstrip.  In  the  general  case  in  which  the  microstrip 
lines  on  either  side  of  the  gap  are  of  different  widths,  this  analysis  must  be 
performed  for  each  side.  Far  from  the  discontinuity,  the  current  on  the  strip  is 
assumed  to  consist  only  of  the  fundamental  mode.  The  current  in  the  uniform 
region  on  the  source  side  of  the  discontinuity  is  represented  by  an  incident  and 
a  reflected  wave  with  an  unknown  amplitude,  T. 

juntform( x  z)  _  ^aijz,(x)(exp(_j/322)  _  r  exp(-|-;73*z)) 

(=1 

Jzn,foTm{X,Z)  =  biJxi(x)(exp(—j 0zz)  +  r  exp (+jfizz)) 

t=i 

The  current  in  the  uniform  region  on  the  load  side  of  the  discontinuity  is  rep¬ 
resented  by  a  transmitted  wave  with  an  unknown  amplitude,  T. 

1=1 

1=1 
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gap 


Figure  6.2:  Arrangement  of  Basis  Functions  for  Analysis  of  Gap  Discontinuity 

The  functions  Jzi  and  Jxi  are  the  Maxwellian  basis  functions  defined 
in  Chapters  3  and  5.  The  constants  at . . .  ,  b\ . . .  bL  and  /?*,  the  propagation 

constant  are  determined  in  the  spectral  domain  analysis  for  uniform  line.  The 
longitudinal  dependance  of  the  current  in  the  perturbed  region  near  the  dis¬ 
continuity  is  augmented  with  piecewise  linear  subsectional  basis  functions.  A 
typical  arrangement  is  shown  in  Figure  6.2  The  transverse  dependance  of  the 
current  in  the  perturbed  region  is  represented  with  Maxwellian  basis  functions 
with  variable  coefficients  thus  taking  into  account  mode  conversion.  Note  that 
the  basis  functions  are  really  those  of  two  microstrip  open-end  discontinuities. 
Therefore,  by  using  the  same  testing  scheme  as  was  used  in  the  analysis  of  the 
open-end  discontinuity  we  can  generate  a  system  of  linear  equations  with  the 
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following  form 


’  [ Kleft ]  [KCTOsa] 

sleft] 

[Jfcroaa ] 

fright  j 

Where  Kle is  the  coefficient  matrix  of  the  linear  system  developed  for  the 
open-end  discontinuity  in  Chapter  4,  KT'aht  is  a  similar  matrix  for  the  right- 
hand  open-end  discontinuity,  KCTOsa  is  a  matrix  developed  from  testing  the  right 
hand  side  with  the  left  and  vice-versa,  and  Sle **  and  Sr'9ht  are  vectors  obtained 
by  testing  the  source  functions  with  the  testing  functions  for  the  left  and  right 
open-end  discontinuities.  These  equations  are  solved  to  obtain  T,  T,  and  the 
current  on  the  conductors. 


6.3  Extension  to  Microstrip  Gap-Coupled  Resonators 

The  method  may  be  extended  to  gap-coupled  resonators  with  the 
inclusion  of  more  subsectional  basis  functions  as  shown  in  Figure  6.3 


6.4  Results 


In  Figure  6.6  and  6.7,  the  magnitudes  of  the  transmission  coefficients 
of  two  gap-coupled  half-wave  resonators  are  shown  calculated  in  the  following 
two  ways.  First  the  scattering  matrix  for  an  isolated  gap  is  calculated  using 
the  spectral  domain  method.  Then  the  scattering  matrix  of  the  resonator  is 
calculated  using  the  scattering  matrix  of  the  isolated  gap  along  with  those  of 
a  section  of  uniform  microstrip  line. 
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Second,  the  spectral  domain  method  is  used  to  model  the  entire  res- 
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Figure  6.3:  Arrangement  of  Basis  Functions  for  Analysis  of  Gap-Coupled  Res¬ 
onator 


|T| 


Figure  6.4:  Magnitude  of  Transmission  Coefficient  of  Gap  Discontinuity 
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Figure  6.5:  Phase  of  Transmission  Coefficient  of  Gap  Discontinuity 


Figure  6.6:  Magnitude  of  Transmission  Coefficient  of  Gap-coupled  Resonator 
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Figure  6.7:  Magnitude  of  transmission  coefficient  of  gap-coupled  resonator 

onator  including  coupling  between  the  gaps  as  shown  in  Figure  6.1.  It  can  be 
seen  that  the  coupling  between  the  two  gaps  has  little  effect  when  the  resonance 
frequency  of  the  resonator  is  low.  However,  for  shorter  resonators  the  coupling 
between  the  gaps  shifts  the  resonance  frequency,  lowers  the  Q,  and  lowers  the 
rejection  above  resonance. 

We  conclude  that  the  spectral  domain  method  can  be  effectively  ap¬ 
plied  to  the  modeling  of  gap-coupled  resonators.  Furthermore,  such  analysis  is 
justified  as  coupling  between  the  two  gaps  of  a  gap-coupled  resonator  can  have 
a  significant  effect  on  the  resonant  frequency,  Q  and  rejection  of  the  resonator. 


Chapter  7 


CONCLUSIONS 


7.1  Summary 

A  method  for  the  analysis  of  planar  circuits  on  open  substrates  has 
been  presented.  Through  the  use  of  the  exact  spectral  domain  Green’s  function, 
the  method  includes  the  effects  of  space  wave  and  surface  wave  radiation.  The 
method  has  been  shown  to  be  applicable  to  1-port  and  2-port  structures.  In 
particular,  the  method  has  been  applied  to  microstrip  open-end  discontinuities 
and  slot-line  short  circuit  discontinuities.  The  method  has  been  applied  to  gap 
discontinuities  and  gap-coupled  resonators. 

7.2  Suggestions  for  Further  Applications 

7.2.1  Extension  to  Alport  Networks 

The  deterministic  spectral  domain  method  presented  here  is  very  gen¬ 
eral  and  should  be  useful  for  analyzing  multiport  planar  networks  on  open 
substrates.  As  was  pointed  out  in  Chapter  6,  the  travelling  wave  formulation 
is  more  efficient  for  analyzing  two- port  networks  than  other  methods  such  as 
[19].  Furthermore,  these  other  methods  are  difficult  to  extend  to  iV-port  net¬ 
works.  To  extend  this  method  to  an  N-port  network  we  simply  excite  one  port 
with  a  known  incident  wave  and  then  formulate  N  unknown  outward  travelling 
waves  as  shown  in  Figure  7.1  for  a  4-port.  Therefore,  in  order  to  determine  the 
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Figure  7.1:  Application  to  N-port  Network 

entire  N  x  N  scattering  matrix  of  a  general  N  port,  N  calculations  must  be 
performed.  Of  course,  symmetry  in  the  network  would  reduce  the  number  of 
necessary  calculations. 


7.2.2  Extension  to  Multilayer  Structures 

By  making  use  of  the  immitance  approach  for  the  Green’s  function 
derivation,  this  method  is  easily  extended  to  multilayer  substrate/superstrate 
structures.  It  would  be  useful  to  study  how  surface-wave  and  space-wave  ex¬ 
citation  can  be  controlled  through  the  use  of  multi-layered  substrates  and/or 
superstrates.  It  should  also  be  noted  that  this  method  may  be  easily  adapted 
through  a  simple  modification  of  the  Green’s  function  to  structures  which  are 
covered  but  laterally  open. 


7.3  Suggestions  for  Improvement 

7.3.1  Asymptotic  Evaluation  of  the  Inner  Products 
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The  most  time  consuming  part  of  the  method  is  the  evaluation  of  the 
asymptotic  portion  of  the  inner  products.  It  is  straightforward  to  determine 
an  asymptotic  form  of  the  Green’s  function  for  large  p  because  the  Green’s 
function  is  separable  in  p  and  9.  However,  since  the  basis  functions  used  here 
are  separable  in  a  and  /?,  an  asymptotic  form  for  large  p  is  not  immediately 
evident.  An  asymptotic  form  for  large  p  would  be  most  useful  because  the 
space  domain  calculation  of  the  asymptotic  portions  of  the  matrix  components 
is  not  compatible  with  the  scheme  given  in  Chapter  3  for  precomputing  and 
storing  the  Green’s  function.  Also  the  space  domain  integration  requires  high 
numerical  precision  due  the  subtraction  of  the  source  and  image  contributions. 
This  problem  is  greater  for  electrically  thin  substrates.  Using  more  rapidly 
converging  basis  functions  does  not  appear  to  be  a  better  approach,  because 
the  basis  functions  must  be  chosen  to  model  the  physical  current  (or  field) 
which  will  necessarily  exhibit  edge  conditions.  These  edge  conditions  in  the 
space  domain  translate  into  high  spatial  frequency  components  in  the  spectral 
domain. 

In  view  of  the  above  difficulties  with  evaluating  the  asymptotic  por¬ 
tion  of  the  integrands,  it  seems  promising  to  solve  the  mixed  potential  integral 
equation  (MPIE)  in  the  spectral  domain  instead  of  the  electric  field  integral 
equation.  This  is  because  the  Green’s  functions  used  in  the  MPIE  converge 
faster  than  those  for  EFIE.  This  is  easy  to  see  by  examining  the  case  of  a 
homogeneous  medium.  The  electric  field  Green’s  function  has  a  singularity 
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at  the  source  point.  It  is  this  source  singularity  that  gives  the  spectral  do¬ 
main  electric  field  Green’s  function  its  asymptotic  tail.  On  the  other  hand  the 
Green’s  functions  used  in  the  MPIE  have  ^  type  singularities  at  the  source 
point.  Therefore,  they  have  more  rapidly  converging  spectral  domain  coun¬ 
terparts.  Although  the  MPIE  requires  the  determination  of  both  charge  and 
current,  the  added  computation  may  be  offset  by  the  savings  gained  by  faster 
convergance.  One  disadvantage  to  this  method  is  that,  as  far  we  have  found,  no 
straightforward  method  such  as  the  immittance  approach  exists  for  determining 
the  necessary  Green’s  functions. 

7.3.2  Generalization  of  Basis  Functions 

Although  in  all  of  the  work  presented  here,  entire  domain  basis  func¬ 
tions  were  used  to  represent  the  transverse  dependence  of  the  current  or  field 
quantities,  it  would  be  useful  to  also  make  use  of  subsectional  basis  functions 
for  both  the  transverse  and  longitudinal  dependance.  For  uniform  line,  en¬ 
tire  domain  functions  such  as  the  Maxwellian  functions  are  most  efficient  for 
modelling  the  transverse  dependance,  just  as  entire  domain  functions  are  most 
efficient  to  describe  the  longitudinal  dependence  of  travelling  wave  functions. 
However,  it  appears  that  in  the  regions  near  discontinuties,  subsectional  ba¬ 
sis  functions  are  more  efficient  for  the  representation  of  the  current  or  field. 
For  example,  in  order  to  analyze  a  step-width  discontinuity,  it  would  be  most 
efficient  to  use  subsectional  basis  functions  to  represent  the  transverse  depen¬ 
dance  of  the  current  near  the  discontinuity.  Also,  it  would  be  useful  to  develop 
more  general  subsectional  basis  functions  for  the  modelling  of  complicated  dis¬ 
continuities  such  as  tapers  and  smooth  bends.  The  use  of  rectangular  basis 
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functions  to  approximate  a  smooth  curve  may  lead  to  error  due  to  the  large 
current  density  along  the  edge.  A  triangular  based  function  would  allow  a  more 
continuous  approximation  of  smooth  curves. 

7.3.3  Computational  Algorithm 

By  far  the  largest  speedups  have  been  obtained  through  modifications 
to  the  computational  algorithm.  It  appears  that  great  improvements  could  still 
be  made  in  this  area.  The  bulk  of  the  computational  effort  occurs  in  evaluating 
the  Green’s  function.  As  was  pointed  out  in  Chapter  3,  a  large  savings  was  real¬ 
ized  by  precomputing  and  storing  the  Green’s  function  so  that  the  entire  matrix 
could  be  integrated  at  once.  This  savings  may  be  increased  by  precomputing 
the  Green’s  function  for  all  /3,  storing  the  values  in  a  database,  and  then  using 
an  interpolation  scheme  to  determine  intermediate  points.  It  appears  that  com¬ 
putational  methods  such  as  the  one  presented  here  are  quite  amenable  to  vector 
processing  because  of  the  large  amount  of  non-interelated  computation.  There¬ 
fore,  despite  the  fact  that  the  method  requires  large  amounts  of  computation, 
it  could  be  utilized  with  small  vector  processors  for  engineering  purposes. 


Appendix  A 

DERIVATION  OF  THE  GREEN’S  FUNCTION  FOR 
THE  GROUNDED  SLAB 


In  this  appendix  we  give  some  of  the  details  of  the  immitance  deriva¬ 
tion  of  the  spectral  domain  Green’s  function.  The  immittance  formalism  is  due 
to  [13]  and  is  included  here  for  completeness.  We  begin  by  noting  that  from 
the  Fourier  transform  relationship 

+oo  +oo 

E(x,y,z)=  —  J  J  E(a,y,fi)exp{-j(ax  +  /3z))dad/3  (A.l) 

-oo  -oo 


the  space  domain  field  and  current  quantities  are  superpositions  of  the  inho¬ 
mogeneous  spectral  waves  propagating  in  the  direction  9  from  the  2-axis  where 
9  =  tan-1  [  ;  r  ).  If  we  define  the  coordinate  transformation 

V  y/a3+0*  J 


u  =  z  sin(0)  —  x  cos(0)  (A. 2) 

v  =  z  cos($)  -f-  x  sin($)  (A. 3) 


we  can  decompose  the  spectral  fields  into  T M-to-y  ( Ey ,  Ev,  and  Hu)  and  T E- 
to -y  ( Hy ,  Hv,  and  Eu).  The  T M-to-y  fields  are  created  by  Jv  and  the  TE-to-y 
fields  are  created  by  Ju.  We  define  transmission  line  characteristic  admittances 
Ytmi  and  Ytm2  as  the  wave  admittances  seen  by  a  transverse-magnetic-to-y 
wave  travelling  in  media  1  and  2  of  Figure  2.1 


Ytmi  = 


-Hi 

Et 


±Hl 

k 


(A.4) 
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YtM2 


-Hj  _  +H~ 

Ei  E- 


(A.5) 


where  the  superscript  implies  waves  travelling  in  the  positive  y  direction 
and  implies  waves  travelling  in  the  negative  y  direction.  Likewise,  we 
define  transmission  line  characteristic  admittances  Ytmi  and  YtM2  as  the  wave 
admittances  seen  by  a  transverse- magnetic- to-y  wave  travelling  in  media  1  and 
2  of  Figure  2.1 


Ytei 

Yte2 


+Hj  _  -H; 
Ei  E- 
+Hj  _  - H ~ 
Ei  E - 


(A. 6) 
(A.7) 


The  y-propagation  constants  may  be  obtained  from  the  dispersion  relations  for 
each  region  and  are  given  by 


7i  =  \J a2  +  P2  -  kl, 

(A.8) 

72  =  \J a2  +  P2  -  tRkl. 

(A. 9) 

The 

wave  impedances  may  be  obtained  from  Maxwell’s 

equations  and  are  given 

by 

v  iueo 

Ytmi  =  - , 

7i 

(A.10) 

-  jue0tR 

Y  TM2  =  i 

72 

(A.l  1 ) 

Ytei  =  — , 
jun 

(A.12) 

v>  72 

YTE2  =  - - • 

(A- 13) 

We  may  now  draw  the  equivalent  circuits  in  for  the  TE  and  TM  modes  as 
in  Figure  A.l.  We  define  the  driving  point  impedances  Y£  and  F/  as  the 
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Figure  A.l:  Equivalent  Networks  for  TE  and  TM  Modes 


admittances  looking  upward  and  downward  at  y  —  d  respectively. 


Yxe  = 


Ye  — 
1 2  — 


yf  = 


Yh  — 
I2  ~ 


■H 


ui 


Ev\ 

+  ^u2 

Ev  2 
+HV  i 
Eui 
-Hv  2 


u2 


At  the  interface,  y  =  d,  the  magnetic  field  is  discontinuous  because 
electric  current  at  the  interface. 

Hui  H  u2  —  J v 
Hv2  Hv1  =  Ju 


(A. 14) 
(A. 15) 
(A. 16) 
(A. 17) 
of  the 

(A. 18) 
(A. 19) 


J 
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At  the  interface,  y  =  d,  the  electric  field  is  continuous. 

Ev\  =  Ev  2  =  Ev 


Therefore  we  can  write 


Eul  —  EU2  —  Eu 


Ev  = 

Y{  +  Y2' 
=  -JvZe 

E  =  --L- 
“  Yf  +  Y2h 

=  -Jjk 


where 


Y{  =  YTm 

Y2  -  yrM2C0th(72d) 

Yxh  =  YTEi 

Y2h  =  YTE2coth(l2d). 

The  transformation  back  into  the  x,  z  coordinates  is 

Ez  =  -  (N2Ze  +  N2xZh)  Jz 
-NZNX  (Ze  -  Zh )  Jx 
Ex  =  —NZNX  ( Ze  -  Zh )  jz 
-  (N2xZe  +  N2Zh)  Jx 


where 


N*  “ 

Nz  =  -JL= 


(A. 20) 
(A.21) 


(A. 22) 


(A. 23) 


(A.24) 
(A.25) 
(A.26) 
(A. 27) 


(A.28) 


(A.29) 


(A. 30) 


(A.31) 
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Therefore,  spectral  domain  Green’s  function  is  given  by 


Zzz  =  ~{N]r  +  N\Zh ) 

(A. 32) 

Zzx  =  -NzNx{Ze  -  Zh) 

(A. 33) 

Zxz  =  —NzNx(Ze  -  Zh) 

(A.34) 

Zxx  =  -{NlZK  +  N]Ze) 

(A. 35) 

Appendix  B 


DERIVATION  OF  THE  GREEN’S  FUNCTION  FOR 
THE  UNGROUNDED  SLAB 


The  derivation  of  the  Green’s  function  appropriate  for  the  analysis  of 
slotline  is  a  simple  modification  of  that  for  the  Green's  function  in  Appendix 
1.  Since  the  medium  below  the  slab  is  the  same  as  that  above,  we  draw  the 
equivalent  circuits  in  for  the  T E  and  TM  modes  as  in  Figure  B.l. 

Therefore,  the  driving  point  impedances  at  y  =  d  are 

yf  - 


V' 

r2 


Yh  — 

1 2  ~ 

r  = 

zh  = 


Ytmi 

(B.l) 

Yte\ 

(B-2) 

Ytm  2 

(  Ytm2  +  Ytmi  coth(72d)\ 

\Ytmi  +  Ytm2  coth(72d)  / 

(B.3) 

YtE2  | 
1 

'  Yte2  +  Ytei  coth(72d)  \ 

( B  .4 ) 

\Yte\  +  Vr£2coth(72d)/ 

(B.5) 

V?  +  Y{ 

i 

(B.6) 

Yxh  +  Yf 

As  pointed  out  in  chapter  4,  the  inverted  Green’s  function  is  then  given  by 

(B.7) 
(B.8) 
(B.9) 


Yxx  =  -N2xi--N]  4- 

r  Ze  Zh 

Yxz  =  -NxNz(j^-j^) 
K„  =  -NM-L--L) 
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Ill 


(B10) 

As  with  the  Green’s  functions  Zzz  . . .  Zxx  for  the  grounded  dielectric 
slab,  Yxx . . .  Yzz  have  a  branch  point  at  p  =  ko.  Unlike  Zzz  . . .  Zxx,  Yxx  . . .  Yzz 
are  not  well  behaved  at  the  branch  point.  This  can  be  seen  by  defining 


Ye 

yh 


_1_ 

Ze 


Y'  +  Y2e 

J_ 

Zh 


Yxh 


+  Y* 


(B.ll) 


(B.12) 


Now  the  Green’s  function  is  rewritten  as 


Yxx  =  -NlYe-N]Yh  (B.13) 

Yxz  =  —NxNz(Ye  -  Yh)  (B.14) 

Yzx  =  -NxNz{Ye  -  Yh)  (B.15) 

Y„  =  -N2zYe-Nl?h.  (B.16) 


It  can  be  seen  that,  because  Ytmi  is  singular  at  p  =  k0 ,  Ye  is  singular  at  p  =  ko . 
It  should  be  stressed,  however,  that  this  point  constitutes  a  branch  point  and 
does  not  contribute  to  the  integral. 

As  with  the  Green’s  functions  Zzz  . . .  Zxx  for  the  grounded  dielectric 
slab,  Yxx . . .  YmZ  do  not  have  a  branch  point  at  p  =  k0. 
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